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near  the  root  because  of  the  obstruction  due  to  the  hub. I .  17 
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IX 


Summary 


The  nature  of  the  flow  around  a  rotating  blade  in  reverse  flow  is  described,  integrating 
results  from  fixed  and  rotating  wing  experiments.  The  highly  3-dimensional  flow  phe¬ 
nomena  do  not  conform  to  expectations  based  on  2-dimensional  airfoil  aerodynamics. 
The  same  NACA0013  blade  was  used  both  as  a  fixed  wing,  and  then  as  a  rotor  blade 
in  a  low-speed  wind  tunnel.  Fixed-wing  results  from  load  measurements  and  flow 
visualization  showed  that  the  sharp-edge  vortex  (SEV)  is  a  primary  feature  when  the 
blade  is  yawed  either  forward  or  backward  in  reverse  flow.  The  aerodynamic  loads  are 
better  modeled  using  the  Polhamus  leading  edge  suction  analogy,  with  vortex-induced 
lift  and  pitching  moment  making  significant  contributions.  Surface  tuft  visualization 
on  the  blade  operated  as  a  fixed  wing,  strengthened  the  evidence  for  such  vortices, 
including  their  bursting  where  they  encountered  adverse  pressure  gradient  near  the 
blade  root.  Particle  image  velocimetry  was  performed  in  the  wind  tunnel  with  the 
blade  held  at  discrete  settings  of  yaw  and  pitch  corresponding  to  expected  rotor  blade 
positions.  The  presence  of  a  strong  SEV  was  confirmed  by  velocity  field  data.  The 
blade  was  next  operated  as  one  of  2  blades  of  a  teetering  rotor,  and  stereo  Particle 
Image  Velocimetry  data  were  acquired  at  azimuths  240  (forward  yaw),  270  (zero  yaw) 
and  300  (backward  yaw)  degrees.  At  240  degrees,  the  vortex- induced  pressure  gradi¬ 
ent  induces  an  inward  radial  flow  overcoming  centrifugal  effects,  but  away  from  the 
vortex,  outward  radial  flow  was  evident.  A  strongly  three-dimensional  and  attached 
SEV  is  evident  under  the  blade  at  240  degrees  azimuth.  This  detaches  and  convects 
with  the  flow,  remaining  close  to  the  blade  by  270  degrees.  The  vortex  seen  at  300 
degrees  is  clearly  detached,  but  growth  of  the  core  diameter  corresponding  to  vortex 
bursting,  causes  strong  suction  under  the  blade.  The  flow  around  the  blunt  edge 
is  again  strongly  3- dimensional,  so  that  the  expectation  of  2-D  flow  separation  was 
debunked.  Some  evidence  of  intermittent  separation  is  seen,  but  the  azimuth-resolved, 
ensemble-averaged  flow  is  mostly  attached  around  the  blunt  edge.  Preliminary  static 
pressure  contours  were  derived  from  the  measured  3-component  velocity  field  using  a 
simple  Bernoulli  equation  conversion  in  this  incompressible  flow,  with  the  effect  of 


x 


LIST  OF  TABLES 


xi 


unsteady  pitch  and  plunge  being  minor.  These  pressure  data  enabled  interpretation  of 
the  surface  pressure  field,  and  the  relative  roles  of  vortex-induced  suction  and  radial 
stress.  In  summary,  it  is  apparent  that  the  reverse  flow  regime  should  be  modeled  and 
analyzed  as  a  case  of  SEV  formation  under  the  very  sharply  swept  blade  immediately 
after  180  degrees  azimuth.  The  SEV  evolves  as  the  sweep  decreases  with  azimuth 
increasing.  In  the  regime  before  240  degrees,  an  attached,  strengthening  SEV  may 
be  expected.  At  some  moderate  sweep  (azimuth  beyond  240  degrees  in  our  case) 
the  vortex  bursts  and  detaches  from  the  surface.  Thereafter  it  convects  with  the 
blade,  but  induces  strong  pressure  effects  on  the  blade  surface  even  as  far  as  300 
degrees  azimuth.  The  blunt  edge  flow  is  highly  3-dimensional,  and  has  much  less 
flow  separation  and  unsteadiness  than  might  be  predicted  from  2-dimensional  airfoil 
aerodynamics. 
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Nomenclature 

b  span  of  the  blade,  m 

c  chord  of  the  blade,  m 

k  reduced  frequency 

R  blade  tip  radius,  m 

Re  chord  Reynolds  number 

r  radial  location,  m 

Uoo  freestream  velocity,  m/s 

Ue  effective  velocity,  m/s 

Ur  radial  velocity,  m/s 

Ut  blade  tip  speed,  m/s 

u  velocity  in  x  direction,  m/s 
v  velocity  in  y  direction,  m/s 
w  velocity  in  z  direction,  m/s 
x  chordwise,  referred  to  the  blade,  m 
y  vertical,  referred  to  the  blade,  m 
z  radial,  referred  to  the  blade,  m 
T  circulation,  m2 / s 

H  advance  ratio 

rotor  angular  velocity,  rad/s 
u  vorticity,  1/s 

if}  azimuthal  angle,  deg 


Chapter  1 

Introduction  and  Synopsis 


1.1  Project  History 

The  retreating  side  of  a  helicopter  rotor  in  forward  flight  poses  several  performance- 
limiting  phenomena.  The  most  famous  of  these  is  the  dynamic  stall  problem.  In 
order  to  balance  the  rolling  moments  between  the  two  sides  of  the  rotor,  the  blade 
pitch  is  increased  as  the  blade  moves  through  the  retreating  part  of  the  disc.  This 
results  in  substantial  portions  of  the  blade  experiencing  stall.  Where  the  pitch  rate 
is  high,  the  blade  retains  dynamic  lift  well  past  the  static  stall  angle  of  attack,  and 
then  experiences  dynamic  stall,  with  sharp  excursions  in  pitching  moment.  The 
precise  timing  of  the  different  phases  of  dynamic  stall  have  proven  to  be  challenging 
to  predict.  Some  of  the  difficulties  stem  from  the  fact  that  rotor  dynamic  stall  is  far 
from  being  a  2-dimensional  phenomenon.  Hence  experiments  and  computations  with 
non-rotating  pitching,  plunging  blades  miss  some  significant  physical  phenomena  that 
are  of  importance  in  the  rotating  case.  Rotation  causes  several  effects.  One  is  the 
obvious  variation  of  the  onset  velocity  from  root  to  tip.  The  second  is  the  occurrence 
of  strong  radial  flow  along  the  blade  surface. 

This  project  started  with  a  study  of  the  radial  flow  on  rotor  blades  in  dynamic  stall. 
From  basic  expectations,  radial  flow  should  have  very  strong  effects  on  blade  lift 
according  to  the  expectation  of  monotonically  increasing  radial  stress.  In  fact  it  does 
not  have  such  strong  visible  effects.  This  was  the  central  mystery,  and  the  challenges 
in  its  resolution  included  those  of  capturing  the  velocity  field  in  the  dynamic  stall 
flowfield.  We  showed  that  the  radial  jet  layer  immediately  above  the  surface  in 
the  stalled  regime,  breaks  up  into  discrete  vortical  shear  layer  structures  with  some 


1 


CHAPTER  1.  INTRODUCTION  AND  SYNOPSIS 


2 


consistent  frequency.  A  similar  phenomenon  occurs  in  the  canonical  problem  of 
flow  induced  over  a  spinning  disk,  when  the  disk  has  edgewise  flow  at  an  angle  of 
attack  superposed.  From  such  studies,  the  frequency  and  strength  of  the  discrete 
structures  could  in  principle  be  predicted,  and  this  could  in  turn  bring  closure  to  the 
3-D  dynamic  stall  problem.  We  also  performed  stereo  Particle  Image  Velocimetry 
(SPIV)  capture  of  the  3-component,  phase-resolved,  ensemble-averaged  (meaning  the 
periodic  component  of  the)  velocity  field  of  the  2-bladed  NACA0013  teetering  rotor  in 
dynamic  stall,  from  inception  of  stall  through  reattachment,  at  several  radial  stations. 
This  is  an  extensive  and  unprecedented  knowledge  base.  Our  project  on  this  problem 
was  completed  with  the  doctoral  thesis  of  Vrishank  Raghav  |4|  in  2014.  Beyond  that 
work,  the  rotor  used  in  the  dynamic  stall  work  was  operated  under  conditions  where 
reverse  flow  occurred.  In  this  problem  the  interesting  phenomena  are  under  the  blade. 
This  report  focuses  on  the  reverse  flow  problem.  This  chapter  is  used  both  as  an 
Introduction  and  as  a  Synopsis,  summarizing  the  significant  results  from  both  the 
dynamic  stall  radial  flow,  and  the  reverse  flow  studies.  The  rest  of  this  report  goes 
into  detail  on  the  reverse  flow  experiments. 


1.2  Introduction  to  Reverse  Flow 


Where  the  advance  ratio  (/i,  the  ratio  of  flight  speed  to  rotor  tip  speed)  begins  to 
approach  unity,  reversed  flow  occurs  on  large  portions  of  the  retreating  blade  of  a 
rotor  disc  in  edgewise  flight,  as  shown  in  Figure  0  This  can  occur  on  helicopter 
rotors  at  high  flight  speed,  the  slowed  rotor  of  an  autogiro  or  compound  helicopter, 
and  on  a  wind  turbine  encountering  a  gust  from  the  side.  The  aerodynamics  of  the 
blade  in  the  reversed-flow  regime  pose  several  uncertainties.  For  instance,  the  center 
of  pressure  shifts  suddenly  from  around  the  25%  to  75%  of  the  chord  as  measured 
from  the  leading  edge  on  the  inboard  parts  of  the  blade  span  that  first  experience 
reverse  flow,  while  remaining  near  0.25%  on  the  outboard  parts.  The  lift  direction 
switches  from  up  to  down  on  the  reversed  flow  part.  These  induce  both  a  quasi-steady 
change  in  bending  and  twisting  moments  on  the  blade,  and  vibratory  loads  as  this 
happens  through  part  of  every  rotor  revolution.  These  loads  may  drive  rotors  into 
instability.  High  and  fluctuating  drag  is  expected  because  of  the  separation  around 
the  blunt  edge  of  the  blade.  With  the  advent  of  rotorcraft  employing  counter-rotating, 
co-axial  rotors,  operation  in  the  reverse  flow  regime  is  becoming  increasingly  common. 


This  report  integrates  results  from  an  experimental  study  of  the  flowheld  and  aerody¬ 
namics  of  rotor  blades  in  reverse  flow.  As  summarized  in  Figure  1.2  prior  approaches 
to  reverse  flow  were  from  a  two-dimensional  airfoil  /  high  aspect  ratio  wing  viewpoint, 
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Figure  1.1:  Variation  of  the  extent  of  the  reverse  flow  region  on  the  retreating  blade, 
with  advance  ratio  // 
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applying  yaw  and  aspect  ratio  corrections  to  2D  airfoil  data.  Such  models  predict 
that  the  region  below  the  blade  features  a  dead  2-D  recirculation  zone,  and  the  blunt 
base  has  massive  separation  with  vortex  shedding.  Even  when  sweep  corrections 
are  applied,  history  effects  are  not  captured.  We  argued  that  the  formation  and 
evolution  of  a  strong  3-dimensional  (helical)  vortex  early  on  the  retreating  blade 
side,  would  be  key  to  the  entire  problem.  At  azimuths  well  before  240  degrees,  the 
sharp  edge  of  the  blade  resembles  the  edge  of  a  highly  forward-swept  wing,  causing 
a  strongly  helical  vortex  to  form,  with  significant  axial  flow  in  its  core.  As  azimuth 
increases  the  yaw  decreases,  and  the  subsequent  evolution  of  the  vortex  bears  much 
similarity  to  the  observed  phenomena  on  delta  wings.  We  explored  this  hypothesis 
that  the  genesis  of  the  aerodynamic  loads  was  best  viewed  through  the  perspective  of 
vortex  flows  occurring  on  sharp-edged  swept  wings  at  angle  of  attack.  A  rectangular 
untwisted  NACA0013  rotor  blade  of  (semi-span)  aspect  ratio  3.47  was  used  as  a 
fixed  wing  in  a  wind  tunnel,  operated  through  a  large  variation  in  yaw  and  smaller 
variations  in  angle  of  attack,  in  both  the  forward  and  reverse-facing  orientations. 
With  the  root  cutout  included  during  operation  as  a  rotor  blade,  this  corresponds 
to  an  aspect  ratio  beyond  8  in  airplane  wing  terms.  The  phenomena  in  reverse  flow 
are  dominated  by  the  sharp-edge  vortex  which  is  generated  along  the  sharp  edge 
of  the  blade.  This  vortex  is  three-dimensional  with  high  values  of  axial  velocity,  ft 
originates  at  the  tip  in  the  case  of  the  forward  quadrant,  and  at  the  root  in  the  aft 
quadrant  of  the  rotor.  This  flowfield  does  not  conform  to  2-D  airfoil  aerodynamics, 
nor  to  the  high-aspect-ratio  ’corrections’  derived  from  Prandlt’s  lifting  line  theory 
approximations,  even  for  a  full-scale  rotor  blade.  The  experiment  here  aims  to  reveal 
the  detailed  flow  phenomena,  as  opposed  to  producing  scale-model  results.  The  rotor 
radius  of  2.13m  was  limited  by  the  2.74m  width  and  2.13m  height  of  the  wind  tunnel 
where  the  rotor  is  operated.  The  larger  chord  enables  the  Reynolds  number  based  on 
chord  to  go  well  above  the  laminar  regime.  It  enables  better  resolution  in  capturing 
the  velocity  and  eventually  the  pressure  distributions  on  the  blade. 


1.3  Prior  Work  on  Reverse  Flow 

The  advance  ratios  experienced  by  the  rotor  blades  of  modern  helicopters  such  as  the 
Sikorsky  X2  and  Eurocopter  A"3  can  exceed  0.8  (5j[6|.  Here,  advance  ratio  is  defined 
as  the  ratio  of  forward  flight  speed  to  the  rotational  tip  speed  of  the  rotor.  At  these 
high  advance  ratios,  significant  sections  of  the  retreating  blade  experience  reverse  flow, 
characterized  by  the  freestream  hitting  the  geometric  trailing  edge  of  the  blade  first 
and  traveling  toward  the  geometric  leading  edge.  The  resultant  flow  field  is  not  well 
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Figure  1.2:  Hypotheses  considered  for  explaining  reverse  flow  aerodynamics:  the 
bottom  one  is  onr  vortex-flow  hypothesis 
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understood  and  is  characterized  by  early  flow  separation,  negative  lift,  and  periodic 
vortex  shedding  |7|.  Reverse  flow  is  a  major  limiter  in  the  design  of  high-speed 
rotorcraft.  Helicopter  designs  that  attempt  to  mitigate  the  effects  of  reverse  flow  date 
back  to  the  early  1970’s.  Fairchild  Republic  Division’s  |8j  Reverse  Velocity  Rotor 
concept  featured  a  negative  pitch  angle  on  the  retreating  blade  in  order  to  achieve 
a  positive  angle  of  attack  in  the  reverse  flow  regime  and  airfoils  with  blunt  leading 
and  trailing  edges.  The  more  recent  Sikorsky  X2  Technology  Demonstrator  |9|  has 
coaxial  rotors  based  on  the  Advancing  Blade  Concept  10  platform  and  also  features 
blunt  trailing  edge  airfoil  sections  near  the  blade  root. 

As  rotor  advance  ratio  rises,  the  region  of  reverse  flow  on  the  rotor  disk  expands, 
reaching  the  blade  tip  at  270°  azimuth  at  an  advance  ratio  p  ~  1.  Prior  work  has  been 
surveyed  in  |2  and  is  summarized  here  by  permission  of  the  authors  (we).  The  most 
basic  method  was  to  take  airfoil  data  (7,11-14  from  2D  wind  tunnel  tests  in  the  region 
around  180  degrees  angle  of  attack,  and  correct  for  yaw  with  a  sinusoidal  correction, 
and  if  needed  an  aspect  ratio  correction.  At  the  next  level,  with  experiments  revealing 
evidence  of  vortices  in  the  flowfield,  ’vortex  shedding’  was  postulated;  both  from 
the  blunt  edge  when  using  NACA0012  type  airfoils  with  blunt  leading  edges,  and 
from  the  sharp  edge  as  from  a  thin  flat  plate  held  at  negative  angle  of  attack  |7|. 
Karman  vortex  street  analogues  were  also  developed.  These  models  have  recently 
been  adapted  to  include  unsteady  pitching  effects,  again  from  airfoil  experimental 
data.  Other  studies  with  unswept  airfoils  have  modeled  the  sharp-edge  flow  as  a 
diffused,  separated  shear  layer,  with  the  blunt  edge  shedding  vortices.  Others  have 
postulated  a  ‘reverse-chord  dynamic  stall’  process  |6[,  explaining  the  occurrence  of 
sharp  pitching  moments  as  well  as  strong  vortices  seen  in  the  flowfield  in  experiments. 


Wheatley  et  al.  1 15 1  performed  force  measurements  on  a  Pitcairn  PCA-2  autogyro 
rotor  at  various  pitch  settings  and  advance  ratios  up  to  /i  =  0.7  and  found  a  negative 
correlation  between  lift  coefficient  and  advance  ratio.  Charles  et  al.  1 16 1  tested 
a  UH-1D  rotor  at  advance  ratios  up  to  //  =  1.1  and  found  that  rotor  performance 
predictions  broke  down  at  p  >  0.5.  They  also  observed  flapping  instability  and  a 
long  transient  response  to  control  input  at  p  =  1.1.  MacCloud  et  al.  117  performed 


tests  on  a  teetering  rotor  with  a  NACA  0012  airfoil  at  advance  ratios  up  to  1.0  and 


observed  a  drop  in  lift  coefficient  at  high  advance  ratios.  Harris  et  al.  18  performed 


correlations  studies  on  these  early  datasets  in  2008  and  concluded  that  OVERFLOW-2 
could  not  accurately  predict  lift,  drag  and  pitching  moments  for  reverse  airfoils  and 
could  not  be  recommended  for  use  at  an  advance  ratio  beyond  0.35.  Recent  work 
includes 


18 


24 


The  survey  by  Quackenbush  1 19 1 ,  though  intended  for  computations, 
captures  many  of  the  phenomena  to  be  expected  from  the  geometry,  and  even  shows 
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a  sharp-edged  vortex.  Harris  |18|  summarizes  high  advance  ratio  rotorcraft  flight 
experience.  Niemiec  1 23 1  and  Carter  1 24 1  describe  airfoils  intended  for  the  reversed- flow 
regime.  Uncertainties  in  lift,  drag,  pitching  moment,  blade  bending  and  twist,  and 


rotor  stability  are  cited  in  [25,  26 1 .  Slowed  compound  rotorcraft  with  other  lifting 
surfaces  are  reported  in  (8, 15  17, 27-33  ,  and  modern  designs  in  |5j34|. 

A  recent  study  on  a  slowed  UH-60A  rotor  performed  by  Datta  et  al.  |6|  showed 


evidence  of  reverse  chord  dynamic  stall  and  large  pitch-link  load  impulses  on  the 
retreating  side  of  the  rotor  disc  at  advance  ratios  near  unity.  Kottapali  |35|  made 
initial  attempts  to  predict  the  blade  loads  from  this  study  using  CAMRAD  II,  but 
they  showed  serious  issues  at  p  >  0.8.  Predictions  made  by  Yeo  36  using  CAMRAD 

performed 


II  showed  fair  airload  and  structural  load  correlation.  Potsdam  et  al. 
coupled  CFD  and  comprehensive  analysis  on  this  dataset  and  predicted  unconventional 
wake  patterns  and  a  lower  surface  vortex  on  the  retreating  blade,  attributing  that  to 
dynamic  stall.  Pitching  moment  predictions  on  the  advancing  and  retreating  sides 
were  not  encouraging.  Lee  et  al  performed  time-averaged  force  measurements  and 
flow  visualization  on  various  airfoil  sections  in  reverse  flow  and  found  a  drag  jump 
at  a  =  180°  due  to  the  unsteady  formation  and  convection  of  a  large  vortex  in  the 
wake.  Lind  et  al.  |7|  also  conducted  reverse  flow  studies  on  static  sharp  and  blunt 
trailing  edge  airfoils.  They  found  that  the  large,  negative  angles  of  attack  of  the 
inboard  section  of  a  rotor  in  reverse  flow  causes  flow  separation  and  the  onset  of 
vortex  shedding.  It  was  suggested  that  the  use  of  a  blunt-trailing  edge  airfoil  with  a 
relatively  linear  lift  curve  slope  would  be  ideal  for  the  reverse  flow  regime.  However, 
Lind  et  al.  did  not  study  the  reverse  flow  aerodynamics  of  static  yawed  airfoils  or 
finite  wings. 


In  summary,  the  past  approaches  to  the  rotor  reverse  flow  problem,  have  focused,  as 
in  the  case  of  dynamic  stall  decades  earlier,  on  the  2-D  airfoil  model.  The  270-degree 
point  is  taken  as  the  starting  point  for  analyses,  where  the  blade  has  no  yaw.  To 
explain  the  vortical  structures  reported  by  experimenters,  computational  researchers 
have  focused  on  pitching  and  plunging  airfoils,  even  venturing  into  ’reverse  dynamic 
stall’.  The  presumption  appears  to  have  been  that  the  rest  of  the  retreating  blade 
sector  can  be  analyzed  as  perturbations  of  the  270-degree  case,  with  yaw  corrections 
applied. 
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1.4  Vortex  Flow  Model  for  Reverse  Blade  Aerody¬ 
namics 


We  suggested  that  approaches  derived  from  airfoil  aerodynamics  were  not  at  all 
appropriate  for  reverse  flow,  even  less  than  for  dynamic  stall.  Our  view  of  reverse  flow 
starts  when  the  rotor  blade  is  just  past  180  degrees  azimuth.  The  sharp  trailing  edge 
of  the  blade  is  now  like  a  highly-swept  strake  or  wing  leading  edge  at  (negative)  angle 
of  attack.  It  generates  a  strong,  tight  vortex.  The  rest  of  the  retreating  blade  sector 
sees  this  vortex  evolving  in  various  ways,  as  discussed  below.  This  hypothesis  was 
first  presented  by  Michael  Mayo,  then  a  first-year  MS  student,  at  the  first  round  of 
the  Lichten  Competition  in  January  2013  at  Georgia  Tech.  He  encountered  a  hostile 
environment.  In  the  year  following  that,  Michael  Mayo  and  his  team  at  our  laboratory 
pursued  the  vortex  flow  hypothesis.  This  model  has  now  been  clearly  proven,  and 
explains  all  observed  phenomena.  In  contrast,  models  based  on  2-D  aerodynamics 
fail  to  predict  even  the  situation  at  270  degrees. 


Our  hypothesis  of  sharp-edge  vortex  flow  is  shown  at  the  bottom  of  Figure  1.2  We 


postulated  that  the  flow  around  and  under  a  sharp,  yawed  edge  behaves  similarly  to 
that  around  and  over  a  sharp-edge  delta  wing  or  other  highly  swept  wing  at  angle 
of  attack.  A  strong,  coherent,  helical  sharp-edged  vortex  is  formed.  Where  this 
encounters  an  adverse  pressure  gradient,  as  near  a  hub,  it  can  ’burst’  and  lift  away 
from  the  surface.  The  vortex-induced  lift  and  pitching  moment  can  be  computed 
using  Polhamus’  suction  analogy 
at  higher  angles  of  attack. 


37,38  ,  as  well  as  nonlinear  models  for  vortex  flows 


A  digression  into  3-dimensional  aerodynamics  is  appropriate  here.  The  aerodynamics 
of  a  general  lifting  surface  [39 1  can  be  represented  in  a  potential  flow  analysis. 
The  potential  flowfield  is  a  solution  to  the  Laplace  equation,  a  statement  of  mass 
conservation.  The  velocity  field  at  any  point  is  the  superposition  of  the  velocities 
’induced’  by  each  of  the  ’singularities’  used  to  represent  the  lifting  surface,  such 
as  sources,  sinks,  doublets  or  vortex  loops  whose  individual  flowfields  are  each  a 
solution  of  the  Laplace  equation.  In  unsteady  flow,  vorticity  in  the  wake  is  also 
represented.  The  strength  of  each  "singularity"  is  found  from  the  condition  of  zero 
normal  velocity  or  the  kinematic  velocity  of  the  surface.  The  solution  of  the  set  of 
equations  representing  the  surface  points  has  two  important  regimes  of  simplification: 
high  aspect  ratio  (Prandtl’s  Lifting  Line)  and  very  low  aspect  ratio  (slender  wing 
theory).  Prandtl’s  lifting  line  approach  realizes  that  gradients  of  flow  properties 
along  the  chord  of  a  high-aspect  ratio  wing  are  much  greater  than  those  along  the 
span.  This  is  used  to  approximate  chordwise  variation  by  a  vortex.  The  spanwise 
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variation  of  strength  of  this  ’bound  vortex’  is  accompanied  by  streaming  vorticity 
whose  ’strength’  is  the  negative  of  the  gradient  of  the  vortex  strength  of  the  lifting 
line.  This  leads  to  a  model  of  a  high-aspect  ratio  wing  as  a  lifting  line  (bound  vortex), 
strong  tip  vortices,  and  an  inboard  vortex  sheet.  The  bound  vorticity  strength  at  any 
spanwise  section  is  related  to  2D  airfoil  aerodynamics. 

At  the  other  extreme  is  Slender  Wing  Theory  for  thin  triangular  wings  of  extremely 
small  aspect  ratio.  Spatial  derivatives  along  the  freestream  direction  are  much  less 
than  those  across  the  freestream.  The  domain  of  influence  affecting  properties  at  a 
given  cross-section,  is  limited  to  local  conditions  and  those  upstream:  downstream 
influence  is  negligible.  The  resulting  pressure  distribution  has  a  sharp  negative 
peak  along  the  high-swept  leading  edges  of  the  wing.  R.T.  Jones  (40]  developed  an 
analogue  between  the  aerodynamics  of  a  given  cross-section  of  a  slender  wing  at  angle 
of  attack  in  steady  flow,  and  an  airfoil  executing  plunge  motion  in  still  fluid.  The 
’non-circulatory  lift’  experienced  by  the  latter  is  equated  to  the  rate  of  change  of 
momentum  of  the  ’apparent  mass’  of  a  cylinder  of  fluid  of  diameter  equal  to  the  chord, 
and  unit  length  along  the  span.  For  a  triangular  wing  (meaning  swept  leading  edge) 
this  rate  of  change  of  momentum  arises  for  an  uncambered  wing  in  steady  flow,  from 
the  sweep.  With  this  formulation,  the  pressure  distribution  over  the  leeside  surface  of 
a  slender  wing  may  be  calculated  from  potential  flow.  As  angle  of  attack  increases 
and  sweep  decreases,  the  pressure  distribution  deviates  from  the  ideal  slender  wing 
prediction  and  evolves  into  that  for  a  leading  edge  vortex  above  a  swept  edge. 

Polhamus’  suction  analogy  (38]  is  that  the  suction  on  the  leeside  of  a  swept  wing  at 
angle  of  attack  corresponds  to  the  rate  of  change  of  momentum  of  the  flow  which  is 
forced  to  curl  around  the  edge  and  stay  attached.  The  dependence  on  local  properties 
is  consistent  with  the  slender-wing  approximation  to  general  lifting  surface  theory 
starting  with  the  assumption  that  chordwise  derivatives  on  a  slender  wing  would  be 
much  lower  than  spanwise  derivatives.  This  approximation  also  leads  to  R.T.  Jones’ 
slender  wing  theory,  from  where  the  spanwise  load  distribution,  and  the  distribution 
of  pressure  coefficient  can  be  derived.  The  Cp  distribution  has  a  sharp  negative  peak 
near  the  sharp  edge.  At  a,  or  at  lower  sweep  for  the  same  a,  the  peak  moves  away 
from  the  edge  as  the  attached  flow  is  replaced  by  first  an  attached  vortex  and  then  a 
detached  vortex.  Thus  our  recommended  approach  will  transition  from  Polhamus’ 
analogy  to  slender  wing  theory  and  then  to  its  modified  form  where  there  is  first  an 
attached  and  then  a  detached  vortex. 

For  a  delta  wing  at  a  given  angle  of  attack,  as  wing  sweep  decreases,  the  vortex  gets 
weaker  (core  expands).  Burst  and  liftoff  (detachment)  occur  faster.  As  the  retreating 
rotor  blade  moves  from  181  deg.  (sharp  edge  at  the  highest  forward  sweep  of  89  deg. 
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into  the  240  deg.  (sweep  of  50  deg.)  regime,  the  SEV  is  well  formed  and  attached  in 
the  reverse  flow  zone.  As  the  blade  moves  towards  270  deg.,  surely  the  sweep  angle 
becomes  too  low  to  sustain  an  attached  SEV,  so  there  is  detachment,  and  as  we  will 
see,  bursting.  The  burst,  detached  vortex,  however,  convects  with  the  mean  flow,  so 
that  at  high  advance  ratios  near  1,  the  vortex  remains  close  to  the  blade  despite  being 
detached.  The  presence  of  this  vortex  all  the  way  to  300  degrees  appears  to  inhibit 
formation  of  a  new  backward-sweep  vortex,  though  one  may  form  at  azimuths  higher 
than  300  deg.,  after  the  vortex  has  convected  away  from  the  blade. 


1.5  Evolution  of  the  Vortex  Flow  Model 


Our  wind  tunnel  experiments  1 1 , 41 , 42 1  have  proceeded  from  fixed-wing  airloads 
(Figure  1.3)  and  tuft  visualization,  to  Particle  Image  Velocimetry  (PIV)  under  a 


blade  held  fixed  as  part  of  a  rotor  in  a  wind  tunnel  HI,  and  then  to  stereo  PIV  with 


the  rotor  at  high  fi,  as  shown  in  Figure  L4  and  1.5|  These  showed  the  presence  and 
dominant  role  of  the  sharp-edge  vortex.  Mayo  et  al.  |1|  used  the  Polhamus  suction 
analogy  |37|  38)  and  the  measured  lift  curve  slopes  (Figure  [l.6[)  to  predict  lift  on 
the  static  yawed  blade  in  the  reverse  flow  regime  and  found  good  correlation  at  yaw 
angles  representative  of  a  slender  delta  wing.  The  movement  of  the  center  of  pressure 
with  changing  angle  of  attack  and  sweep  in  the  reverse  flow  regime  is  illustrated  in 
Figure  1.7  from  |l].  Tuft  behavior  in  forward  yaw  versus  backward  yaw  showed  stark 
differences  in  flow  topology.  A  sample  is  shown  in  Figure  1.8  from  (l) . 


Figure  1.9  shows  these  SEVs  at  r/R  of  0.5,  at  advance  ratios  0.7,  0.85  and  1.0  (they 
are  all  at  the  same  radial  station),  much  as  on  a  slender  delta  wing.  The  structures 
are  seen  at  different  r/R  as  well,  as  shown  in  our  prior  papers.  Later  experiments 
showed  vortex  ’bursting’  at  inboard  stations.  Figure  1.10  attempts  to  present  a 
simplified  representation  of  what  is  seen  in  Figure  [L9|  In  |44|  we  explored  the  velocity 
held  under  the  rotor  using  SPIV,  and  then  reconstructed  the  static  pressure  on  the 
blade  using  the  incompressible  Bernoulli  equation.  In  the  following  chapters,  the 
experiments  and  results  are  discussed  in  more  detail. 
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Figure  1.3:  Rotor  blade  used  in  fixed-wing  mode 
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Figure  1.4:  Test  section  as  viewed  from  upstream,  showing  stereo  PIV  system  looking 
at  bottom  of  the  retreating  blade 
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xp  =  0° 


Figure  1.5:  Schematic  of  top- view  of  the  rotor  disk  (not  to  scale),  for  stereo  PIV 
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Figure  1.6:  Lift  curve  slopes  in  forward  and  reverse  sweep  for  a  blade  held  as  a  fixed 
wing.  From  [Tj 
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Figure  1.7:  Movement  of  center  of  pressure  on  a  fixed  wing  in  reverse  flow,  with  yaw 
and  angle  of  attack,  at  chord  Reynolds  number  of  300K.  From  JTJ 
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a  =  195‘ 


Figure  1.8:  Surface  tuft  behavior  shows  clear  evidence  of  a  vortex  and  vortex  bursting 
on  a  wing  in  reverse  flow  at  45  degree  forward  yaw  and  15  degrees  angle  of  attack. 
From  [lj 


r/R  =  0.7 


r/R  =  0.4 


Blade  Tip 


Trailing-edge 


Figure  1.9:  Sharp-edge  vortex  captured  using  stereo  PIV  at  r/R  of  0.5  under  a 
rotating  blade  at  240  deg.  azimuth,  for  several  advance  ratios.  From  |2 
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Figure  1.10:  A  simplified  artist’s  conception  of  the  SEV  at  240  degrees,  as  it  might  be 
for  a  blade  of  larger  aspect  ratio.  Note  that  the  vortex  starts  at  about  r/R  of  0.7,  and 
grows  inward.  The  dynamic  pressure  of  the  relative  flow  is  highest  near  the  root,  but 
the  vortex  may  also  burst  near  the  root  because  of  the  obstruction  due  to  the  hub. 


Chapter  2 

Experimental  Methods 


2.1  Wind  Tunnel 


Experiments  were  performed  in  the  2.13mx2.74m  John  Harper  wind  tunnel  at  our 
institution.  The  untwisted  rectangular  blade  has  a  NACA  0013  planform.  The  rotor 
used  in  the  rotating  experiments  has  two  such  blades  attached  to  a  teetering  rotor 
hub  actuated  by  a  3.73  kW  motor.  Figure  ??  shows  the  fixed  wing  and  rotor  setups 
with  the  laser  sheet  used  for  stereo  PIV.  Please  refer  to  (3  43  for  further  detail. 


The  same  blade  was  used  in  fixed-wing  experiments  over  a  wide  range  of  yaw  and  angle 
of  attack  in  both  forward-facing  and  reverse  flow,  with  6-DOF  load  measurements 
and  tuft  visualizations  conducted.  The  blade  was  then  mounted  on  the  rotor,  and 
held  stationary  at  discrete  values  of  yaw  corresponding  to  rotor  azimuths  of  240, 
270  and  300  degrees,  at  angles  of  attack  in  the  range  of  those  encountered  in  the 
rotor  tests.  Here  it  is  noted  that  the  angle  of  attack  calculation  is  approximate.  In 
actual  rotating  mode,  the  local  section  angle  of  attack  will  vary.  This  is  because 
the  inflow  will  vary  due  to  the  wake  influence,  and  the  local  relative  velocity  due  to 
rotation.  Only  an  expected  average  value  can  be  set  in  fixed-wing  mode.  Particle 
image  velocimetry  was  used  to  study  the  flowheld  under  the  blade  in  reverse  flow. 
These  are  compared  with  results  obtained  later  under  the  rotating  blade. 
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2.2  Instrumentation  and  data  processing 


A  NACA  0013  untwisted  rectangular  blade  with  an  aspect  ratio  AR  =  3.49,  a  chord 
c  =  0.178  m,  and  a  span  b  =  0.622  m  was  set  up  as  shown  in  Fig.  m  Stepper- 
motors  were  used  to  control  pitch  and  yaw.  A  6-DOF  load  cell  was  used  to  measure 
aerodynamic  loads.  A  5Hz  low  pass  digital  filter  removed  noise,  and  a  calibration 
matrix  converted  the  filtered  voltage  values  into  forces  and  moments.  Raghav  et 
al.  |3|  further  details  the  instrumentation  and  data  processing.  This  blade  is  part 
of  the  2-bladed  teetering  rotor  system  already  used  in  forward  flight  experiments  in 
this  tunnel,  and  also  used  in  experiments  at  high  advance  ratios.  Thus  the  aspect 
ratio  and  root  structure  of  the  blade  were  dictated  by  the  requirements  of  the  rotor 
experiments,  and  it  was  essential  to  use  the  same  in  these  fixed-wing  experiments  for 
a  valid  comparison  with  rotating-blade  aerodynamics.  A  high  torque  stepper  motor 
with  a  resolution  of  0.1°  was  used  to  control  angle  of  attack,  a,  and  a  high  torque 
stepper  motor  with  a  resolution  of  1°  was  used  to  control  yaw  angle,  T.  The  motors 
were  controlled  by  the  uDrive  NI  motion  controller.  A  picture  of  the  rotor  blade 


setup  inside  the  wind  tunnel  is  shown  in  Fig.  2.1 


In  rotating  mode,  the  collective  pitch  angle  was  7°  and  the  longitudinal  cyclic  angle  was 
—8°.  This  results  in  a  15°  pitch  at  T  =  270°.  The  phase  locked  SPIV  measurements  on 
the  rotating  rotor  blade  were  acquired  at  advance  ratios  of  /i  =  0.7, 0.85, 1.0  and  a  rotor 
angular  velocity  Q  =  200rp?77,(20.94rad/s).  Measurements  were  gathered  at  radial 
locations  of  r/R  =  0.4,  0.5,  0.514,  0.6,  0.7  and  at  azimuthal  angles  4/  =  240°,  270°,  300° 


summarizes  the  flow  and  test  conditions. 


equivalent  to  ^static  —  30°,  0°,  —30°  for  the  static  blade  in  prior  work  1 1  42  .  Table  3.4 


Table  2.1:  High  advance  ratio  facility  rotor  specifications 


Parameter 

Value 

Blade  span,  chord  (m) 

0.622,0.178 

Aspect  ratio,  solidity 

3.49,0.0895 

Radius,  floor  height  (m) 

0.889,1.3 

Precone,  max  TPP  tilt° 

1.6,16 

Maximum  collective,  cyclic  ° 

10,6.5 
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Figure  2.1:  Side  view  of  the  rotor  setup  in  the  wind  tunnel,  where  Xb  and  are  blade 
fixed  coordinates) 

2.3  Stereoscopic-PIV  instrumentation 

Stereoscopic-PIV  (SPIV)  was  used  to  investigate  the  3D  flow  field  under  the  rotating 
blade.  A  double-cavity  Nd:YAG  laser  delivered  200  mJ,  5ns  pulses  at  up  to  14 
double-pulses  per  second.  One  velocity  field  was  accessible  per  rotor  revolution, 
synchronized  to  the  same  blade.  The  light  sheet  was  2  mm  thick.  Seeding  with 
approximately  10 /im  droplets  came  from  a  Laskin-nozzle  aerosol  generator.  The 
blade  was  coated  with  matte  black  paint  to  minimize  laser  reflections  from  the  lower 
surface.  Surface  scattering  was  masked  during  velocity  vector  computations.  A  series 
of  100  phase-locked  image  pairs  (implying  over  100  revolutions,  or  30  seconds  at 
200  rpm)  was  acquired  at  each  measurement  location  using  two  LaVision  PRO-X 
2M  cameras  with  50  mm,  f/8  lenses,  1600x1200  pixel  resolution  and  pixel  size  of 


CHAPTER  2.  EXPERIMENTAL  METHODS 


21 


Table  2.2:  Test  matrix  for  experiments 


Azimuthal  Radial  sta- 

an-  tion  (r/R) 

gle  (4'°) 

Advance  Ra¬ 
tio  (/i) 

240° 

0.4,  0.5,  0.6,  0.68 

0.7,  0.85,  1.0 

270° 

(0.7  for  T  =  240°) 

300° 

7.4x7. 4/im2.  The  stereo  angle  between  camera  optical  axes  varied  between  20°  and 
39°  for  different  radial  locations.  The  particle  size  in  the  camera  image  ranged  from 
1.76  pixels  to  2.81  pixels.  LaVision  DaVis  8.0  software  was  used  to  compute  velocities 
from  spatial  cross-correlation  between  the  images.  An  interrogation  window  overlap  of 
50%  and  a  second  interrogation  pass  with  a  reduced  window  size  were  used  to  increase 
the  signal-to-noise  ratio  of  the  correlation  peak.  The  first  and  second  passes  used 
interrogation  windows  of  64x64  and  32  x32  pixels.  This  yielded  a  spatial  resolution 
between  velocity  vectors  of  1.4  mm  -  2.2  mm  ~  0.008c-0.012c.  An  allowable  range  of 
vector  magnitudes  and  a  median  filter  were  specified  to  reduce  noise  at  the  edges  of 
the  field  of  view.  Figure  |2T2] shows  a  velocity  measurement  window  and  the  orientation 
of  pressure  contour  plots. 


2.4  Accuracy  estimates 


Table  2.3:  Summary  of  velocity  uncertainty  estimates 


Type 

Error 

In-plane  random  error 

0.074-0.252  px 

In-plane  bias  error 

0.01  m/s 

In-plane  measurement  error  ( %UP ) 

0.2%-0.96% 

Out-of-plane  measurement  error  (%WP) 

0.29%-1.57% 

Estimates  are  repeated  from  |4 , 43 1  by  permission,  and  summarized  in  Table  2.3 
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Freestream 

Direction 


Blunt  leading  edge  (BE) 


r/R  =  0.68 


Bottom  Surface  of  Blade 


r/R  =  0.4 


Sharp  trailing  edge  (SE) 

Cp  Contour  plots  location 


-x/c  =  0 


(b)  Bottom  surface  Cp  contour  region 


Figure  2.2:  Measurement  Regions 
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Table  2.4:  Blade  Motion  Amplitude  and  Phase 


Advance 

Amplitude 

A,  Phase, 

Ra¬ 
tio  (/i) 

deg.  (r/R) 

deg  (0) 

0.7 

3.38° 

-31.1° 

0.85 

3.47° 

-20.18' 

1.0 

3.73° 

1 

7U 

o 

2.5  Tip  Path  Plane  (TPP)  and  Blade  Vertical  Mo¬ 
tion 


The  lateral  and  longitudinal  tilt  of  the  TPP  were  measured  using  two  video  cameras 
calibrated  against  a  grid-board.  The  results  were  first  shown  in  [2|.  The  vertical 
velocity  of  the  blade  at  each  station  was  calculated  from  phase-shifted  sinusoidal 
fits  with  the  measured  values  at  0-180  and  90-270.  At  each  /i  and  r/R,  blade  speed 
variation  with  azimuth  is  given  by: 


W  =  A(y^(ylR>nC0S(^ 


0) 


(2.1) 


A  and  0  are  given  in  Table  2.4  The  curve  fits  are  accurate  to  well  under  0.5  percent 
error,  and  the  measurement  error  is  discussed  above,  uncertainty  in  fitting  with  only 
two  points  precludes  converting  velocity  data  to  true  blade-fixed  coordinates.  The 
blade  vertical  speed  in  Figure  [273],  found  by  differentiating  the  curve  fit  for  position, 
can  be  on  the  order  of  lm/s  as  /i  rises.  TPP  orientation  does  not  agree  with  textbook 
speculations.  This  is  not  surprising  because  of  the  uncertain  aerodynamics  in  the 
reverse  flow  region.  In  converting  to  ’blade-fixed’  coordinates  to  clarify  vortical 
structures,  only  the  predetermined  constants  arising  from  free  stream  speed  and  blade 
section  tangential  speed  are  subtracted. 
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Figure  2.3:  Blade  vertical  speeds  in  m/s  for  different  advance  ratios  in  the  reverse 
flow  regime 


Chapter  3 

Fixed  Wing  Results 


In  this  chapter  we  pursue  the  hypothesis  that  the  sharp  trailing  edge  of  a  reverse  flow 
blade  develops  an  attached  vortex  in  a  similar  manner  to  a  sharp-edged  delta  wing. 
Lift,  drag  and  pitching  moment  data  for  a  static  yawed  blade  in  reverse  flow  are 
acquired  using  a  load  cell  and  the  lift  measurements  are  compared  to  predictions  made 
with  an  adaptation  of  the  Polhamus  model.  Surface  tuft  flow  visualization  confirms 
the  presence  of  an  attached  span-wise  vortex  on  the  blade.  Surface  flow  behavior 
is  used  to  explain  the  behavior  of  observed  aerodynamic  loads,  and  correlations  are 
made  to  previous  results  found  in  the  literature.  Vorticity  transport  analysis  is  used 
to  justify  the  surface  flow  behavior  variation  with  yaw  angle  and  angle  of  attack. 

The  first  step  in  this  process  was  the  exploratory  wind  tunnel  study  performed  by 
Raghav  et  al.  |3|  on  reverse  flow  over  a  static  yawed  wing.  They  found  evidence  of 
vortex-induced  lift.  Load  cell  measurements  showed  higher  lift  curve  slopes  for  reverse 
flow  test  cases  when  compared  to  regular  flow  cases  at  the  same  yaw  angles,  and 
tuft  visualization  confirmed  the  existence  of  an  attached  vortex  on  the  rotor  blade  in 
reverse  flow.  In  that  Note  we  presented  a  hypothesis  that  vortex  lift  similar  to  that 
on  the  swept  leading  edge  of  a  delta  wing,  may  be  a  strong  phenomenon  both  in  the 
case  of  the  yawed  reversed-flow  static  blade,  and  in  the  case  of  the  retreating  rotor 
blade. 


In  this  chapter  we  present  results  from  further  analysis  from  these  experiments, 
and  further  experiments  over  a  range  of  yaw  angles.  In  addition,  the  vortex  flow 
hypothesis  from  the  previous  study  is  pursued  and  expanded  by  using  the  Polhamus 
suction  analogy  |37|  for  delta  wing  lift  to  predict  reverse  flow  lift  curves.  Drag  polars, 
pitching  moments  and  chordwise  center  of  pressure  movement  are  also  analyzed  in 
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addition  to  flow  visualization  analysis  at  various  yaw  angles  and  Reynolds  numbers. 
Use  of  a  static  yawed  wing  allows  rotational  effects  to  be  decoupled  from  the  already 
complex  behavior  of  the  reverse  flow  region.  The  wing  used  in  this  study  is  part  of  a 
two-bladed  rotor  used  in  dynamic  stall  research  by  Raghav  1 3 1 .  The  flow  regime  is 
limited  to  that  of  low  Mach  number  and  moderate  chord  Reynolds  numbers. 


3.1  Polhamus  suction  analogy  for  yawed  rotor  blade 


The  Polhamus  suction  analogy 
flow  as  shown  in  Equation 
discussion  section. 


3.1 


is  adapted  for  the  case  of  a  yawed  wing  in  reverse 
The  predictions  are  presented  in  the  results  and 


CL  —  Kp  ■  sina 


cos2  a  +  (  Kv  —  K, 


1  \ 


'  ^ _  I  _ 

p  neAR  J  cosT 


cosa  ■  sin  a 


(3.1) 


The  terms  on  the  right-hand  side  correspond  to  the  lift  generated  by  potential  flow 
and  the  leading-edge  vortex  respectively.  Kp  depends  only  on  the  plan-form  shape 
and  gives  the  lift-curve  slope  from  small-perturbation  theory  for  small  angles  of  attack. 
In  this  paper,  Kp  will  be  set  equal  to  the  lift-curve  slope  of  the  rotor  blade  in  regular 
flow  for  each  yaw  angle  being  examined.  The  delta  wing  sweep  angle  A  is  equivalent 
to  the  yaw  angle  T.  of  the  blade  used  in  this  work.  Oswald’s  efficiency  factor  e  is 


calculated  using  Shevell’s  method  45' . 


3.2  Flow  and  test  conditions 

The  experiments  were  performed  at  two  Reynolds  numbers  of  1.7  x  105  and  3.0  x  10s 
at  yaw  angles  ranging  from  T  =  —60° (backwards  yaw)  to  T  =  60° (forwards  yaw). 
The  test  cases  are  summarized  in  Table  13.41  and  described  below: 

1.  The  first  test  validated  the  experimental  method  |3|  by  comparing  results  for 
regular  flow  against  published  data. 

2.  The  second  test  acquired  lift,  drag  and  pitching  moment  data  for  the  rotor 
blade  in  reverse  flow  at  two  Reynolds  numbers  and  various  yaw  angles.  Angle  of 
attack(a)  resolution  was  lower  in  the  linear  region  of  the  lift  curves  and  higher 
near  the  stall  region. 
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Table  3.1:  Test  matrix  for  experiments 


Test 

Flow  Type 

Data 

Acquired 

T  and  a 

Ae(xl05) 

1 

Regular  Flow 

Validation 

0°  <  T  <  60°,  0°  <  a  <  30° 

1.7 

2 

Reverse  Flow 

Lift,  Drag, 

Pitching 

Moment 

-60°  <  T  <  60°  ,  180°  <  a  <  210° 

1.7,  3.0 

3 

Reverse  Flow 

Flow  Visu¬ 
alization 

T  =  ±30, ±45,  180°  <  a  <  220° 

1.7,  3.0 

3.  The  third  test  consisted  of  surface  flow  visualization  using  tufts  at  two  Reynolds 
numbers  and  various  yaw  angles.  The  rotor  blade  was  pitched  through  the 
angle  of  attack  range  in  increments  of  5°. 

Dynamic  pressure  was  measured  using  a  Barotron  pressure  transducer  connected 
to  a  Pitot-static  probe  in  the  test  section.  Reynolds  numbers  were  set  at  170,000 
and  300,000  by  adjusting  tunnel  speed  as  gauged  from  the  dynamic  pressure,  and 
ambient  temperature  and  pressure  measurements.  Lift  and  drag  were  measured  by 
setting  the  rotor  at  a  desired  T,  the  wind  tunnel  to  a  speed  corresponding  to  the 
desired  Reynolds  number,  and  stepping  the  blade  through  the  desired  angles  of  attack. 
A  T  of  0°  corresponds  to  the  blade  span  being  perpendicular  to  the  freestream.  If 
the  rotor  blade  trailing  edge  is  facing  into  the  flow,  then  T  =  0°  corresponds  to  a 
helicopter  blade  azimuth  angle  of  270°  as  seen  in  Fig.  ??.  Yaw  angles  above  60° 
were  not  used  because  an  advance  ratio  greater  than  2  would  be  necessary  to  achieve 
reverse  flow  in  this  region  of  a  helicopter  rotor  disc.  Reverse  flow  conditions  were 
achieved  by  rotating  the  rotor  blade  180°  so  that  the  trailing  edge  faced  into  the  flow. 
In  this  paper,  a  =  180°  corresponds  to  the  sharp  trailing  edge  of  the  blade  facing  into 
the  flow  with  no  angle  of  attack,  and  a  =  200°  would  correspond  to  the  the  sharp 
trailing  edge  facing  into  the  flow  with  20°  angle  of  attack.  Reverse  flow  data  were 
acquired  in  the  range  of  T  =  —60 °  (backward  sweep )  to  T  =  60°  (forward  sweep ) 
and  at  a  =  180°  —  190°  in  increments  of  2°,  at  a  =  190°  —  196°  in  increments  of  1°, 
at  a  =  198°,  and  at  a  =  200°.  Data  were  acquired  at  a  values  greater  than  200° 
when  necessary  to  observe  the  stall  region.  The  force  and  moment  data  presented  are 
ensemble  averaged  data  with  approximately  300,000  samples  taken  at  each  a. 
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3.3  Error  analysis 


The  total  error  in  measurement  of  the  aerodynamic  load  coefficients  was  computed  to 
be  ±1.61%  in  the  worst  case.  The  total  error  was  computed  using  the  individual  errors 
of  the  parameters  listed  in  Table  ±2  The  ensemble  averaged  data  are  repeatable 
with  a  standard  deviation  within  4.5%  of  the  measured  average  loads.  The  blade 
airfoil  contour  was  measured  at  18  different  locations  across  the  span  of  the  blade  to 
quantify  the  difference  from  NACA0013  coordinates.  The  root  mean  square  error  for 
the  difference  was  less  than  0.1%  of  the  chord. 


Table  3.2:  Summary  of  errors 


Parameter 

Error  Estimate 

Density 

0.03% 

Velocity 

0.04% 

Reynolds  number 

±470 

Load  cell  calibration 

±1.37% 

Yaw  angle 

±0.05° 

Angle  of  attack 

±0.01° 

3.4  Aerodynamic  loads 

3.4.1  Lift 


Lift  curve  slopes  for  the  reverse  flow  regime  at  positive  and  negative  yaw  angle  and 
two  different  Reynolds  numbers  are  shown  in  Fig.  3.1  along  with  forward  yaw  regular 
flow  lift  curve  slopes  at  Re  =  1.7  x  105.  Figure  2.1  defines  the  positive  (forward)  and 
negative  (backward)  yaw  angles.  Lift  curve  slopes  were  calculated  using  the  linear 
region  in  the  equivalent  a  range  of  0°  —  10°.  The  occurence  of  stall  near  a  =  10 


m 


both  regular  and  reverse  flow  regimes  agrees  with  Lind  et  al.  |7|  for  a  NACA  0012 
airfoil  in  regular  and  reverse  flow  and  with  NACA  0012  reverse  flow  data  collected 
by  Critzos  et  al. 


12 


in  the  Langley  low- turbulence  tunnel.  Raghav  et  al.  |3|  showed 
that  the  regular  flow  lift  curve  slopes  follow  the  cos2T  fit  predicted  by  cross-flow 
equivalence 
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The  lift  curve  slopes  for  the  forward  yaw  reverse  flow  cases  are  slightly  higher  than 
the  corresponding  regular  flow  lift  curve  slopes,  which  suggests  that  there  is  an  extra 
source  of  lift  in  the  reverse  flow  regime.  The  difference  in  the  regular  and  reverse 
flow  lift  curves  seen  in  this  study  correlates  well  with  the  difference  in  lift  coefficient 
seen  by  Bartlett  [47]  for  beveled  and  rounded  wings  with  sweepback  angle.  The  lift 
curve  slopes  at  Re  =  3.0  x  105  are  slightly  higher  than  the  ones  at  Re  =  1.7  x  105, 
which  suggests  that  Reynolds  number  does  have  an  effect  on  aerodynamic  loads  for 
a  rotor  blade  in  the  reverse  flow  regime.  The  difference  between  the  two  Reynolds 
numbers  is  small  and  even  non-existent  for  certain  yaw  angles,  which  supports  the 
suggestion  made  by  Lind  et  al.  |7  that  the  aerodynamic  forces  in  reverse  flow  are  less 
dependent  on  Reynolds  number  due  to  early  separation.  Reynolds  number  appears  to 
have  little  if  any  effect  on  the  forward  yaw  reverse  flow  lift  curves.  This  is  consistent 
with  the  negligible  effect  of  Reynolds  number  on  aerodynamic  loads  for  sharp-edged 
delta  wings  |48|. 

The  backward  yaw  lift  curve  slopes  do  not  follow  the  expected  cos2T  fit,  and  Reynolds 
number  does  appear  to  have  more  effect  on  the  lift  curve  slopes  at  T  =  —15°  and 
T  =  —30°  when  compared  to  forward  yaw  angle  cases.  Possible  explanations  for  this 
behavior  are  the  different  orientation  of  the  mount  and  wingtip  in  the  backward  yaw 
case  -  this  is  further  examined  in  the  flow  visualization  section. 


Equivalent  Yaw  Angle  (°) 


Figure  3.1:  Lift  curve  slopes  at  various  T  and  Re  (T  negative  for  backward  yaw) 
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3.4.2  Predictions  using  Polhamus  suction  analogy 


The  Polhamus  suction  analogy  described  in  section  3.1  was  used  to  predict  lift 
curves  for  the  reverse  flow  rotor  blade  at  Re  =  1.7  x  105  and  forward  yaw  values  of 
T  =  0°,  15°,  30°,  45°,  and  60°.  These  yaw  values  were  chosen  because  they  allowed  the 
Kp  term  in  the  vortex  lift  equation  to  be  matched  to  the  lift  curve  slopes  obtained 


from  the  regular  flow  validation  trials.  Figure  3.2(a)  compares  the  predictions  made 
by  the  vortex  lift  model  at  T  =  45°  and  Re  =  1.7  x  105  to  experimental  data.  The 
solid  line  is  the  prediction  made  using  the  full  equation  3.1,  and  the  dashed  line 


is  the  prediction  made  using  equation  3.1  with  Kv  (second  term  on  the  right  hand 
side)  set  to  zero.  The  dashed  line  corresponds  to  lift  predicted  by  potential  theory 
only.  The  full  vortex  lift  model  predicts  the  lift  coefficient  with  good  accuracy  up 
until  a  =  195°,  where  the  experimental  lift  curve  begins  to  more  closely  follow  the 


potential  prediction.  In  reality,  at  higher  angle  of  attack  the  vortex  tends  to  burst  49 


and  the  measured  lift  will  decrease.  However,  the  current  mathematical  model  does 
not  account  for  the  vortex  burst  phenomenon.  This  is  reconciled  in  the  surface  flow 
visualization  section  347  where  vortex  burst  is  observed  at  higher  angle  of  attack.  It 
should  be  noted  that  the  experimental  lift  curve  is  higher  than  the  potential  prediction 
well  after  diverging  from  the  vortex  model,  which  may  signify  that  some  vortex  lift  is 
being  generated  well  after  breakdown  begins. 

Lift  curve  slopes  were  calculated  using  the  lift  curve  predictions  generated  by  the 
vortex  lift  model  and  are  compared  to  experimental  data  in  Figure  3.2(b)  The  vortex 
lift  predictions  are  fairly  accurate  for  4/  >  30°.  The  over- prediction  at  lower  4/  stems 
from  the  fact  that  the  Polhamus  equation  for  low  aspect  ratio  delta  wings  suffers  at 
low  4>  values  which  correspond  to  the  leading  edge  sweep  of  a  high  aspect  ratio  delta 
wing. 


3.4.3  Drag  and  pitching  moment 


Drag  polars  for  the  regular  and  reverse  flow  regime  at  T  =  30°  are  shown  in 
Figure  3.3(a)  The  drag  polars  at  this  4/  are  fairly  representative  of  drag  polars  at 


other  4/.  The  L/D  ratios  are  higher  in  the  regular  flow  cases.  This  agrees  with  results 
for  delta  wings,  where  a  sharp-edged  delta  wing  has  a  lower  L/D  than  a  delta  wing 
with  a  rounded  leading  edge  1 50 1 .  The  regular  flow  L/D  ratio  at  a  =  10°  is  almost 
twice  as  high  as  the  corresponding  reverse  flow  L/D.  This  result  agrees  with  that 
discussed  by  Lind  et  al.  |7]  showing  that  a  NACA  0012  airfoil  in  reverse  flow  would 
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(a)  Lift  curves  for  4'  =  45°  produced  by  vortex  lift(b)  Lift  curve  slopes  obtained  from  vortex  lift 
model  and  experiment  model  and  experiment  at  Re  =  1.7  x  105 


Figure  3.2:  Lift  predictions  using  the  Polhamus  model  for  a  yawed  blade  in  reverse 
flow 


produce  almost  twice  as  much  drag  at  the  stall  angle  while  still  having  the  same 
C 'Lmax  as  a  regular  flow  airfoil. 

The  reverse  flow  pitching  moments  were  proportional  to  lift  in  the  linear  region  of 
the  lift  curve.  A  plot  of  Cl  and  Cm  at  T  =  45  for  Re  =  1.7  x  105  and  3.0  x  105 
is  shown  in  Figure  3.3(b)  The  non-zero,  negative  pitching  moment  in  the  linear 


lift  region  suggests  that  the  center  of  pressure  was  located  closer  to  the  effective 
leading  edge  (i.e.,  the  sharp  edge).  The  chordwise  center  of  pressure  was  computed 
from  the  normal  aerodynamic  force  using  the  lift  and  drag  measurements  as  given  by 


equation  3.2  and  the  measured  pitching  moment. 


N  =  Lcos(a )  +  Dsin(a) 


(3.2) 


It  should  be  noted  that  the  resultant  chordwise  center  of  pressure  location  is  an 
average  over  the  span  of  the  blade.  The  chordwise  center  of  pressure  locations  for 
various  forward  yaw  angles  at  Re  =  1.7  x  105  and  Re  =  3.0  x  105  are  shown  in 
Figures  [3.4(a)|  and  3. 4(b)| respectively.  A  0.75c  location  on  the  plot  corresponds  to 
the  quarter  chord  of  the  reverse  airfoil.  Data  points  derived  from  lift  values  lower 
than  0.3N  were  omitted  due  to  high  uncertainty. 
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Both  plots  show  that  at  low  T,  the  center  of  pressure  is  even  closer  to  the  effective 
leading  edge.  A  possible  reason  for  this  trend  is  discussed  in  the  flow  visualization 
section.  The  approximate  0.75c  center  of  pressure  location  for  both  plots  agrees  with 
data  collected  by  Datta  |6|.  For  Re  =  1.7  x  105,  the  center  of  pressure  locations  start 
near  0.7c  and  move  toward  the  effective  leading  edge  as  a  is  increased.  Beyond  the 
linear  lift  region,  the  center  of  pressure  moves  toward  the  effective  trailing  edge  and 
eventually  settles  near  0.75c  at  a  =  200°.  The  Re  =  3.0  x  105  cases  show  a  different 
trend,  with  center  of  pressure  locations  that  start  closer  to  the  effective  leading  edge 
and  move  toward  the  effective  trailing  edge  as  a  is  increased,  eventually  settling  near 
0.75c.  Reasons  for  the  center  of  pressure  movement  will  be  discussed  in  the  flow 
visualization  section. 


Equivalent  Angle  of  Attack  (°) 


Angle  of  Attack  (°) 


(a)  Drag  polars  at  it  =  30°  (a  =  180°  —  200°  for  (b)  Variation  of  lift  and  pitching  moment  at 
reverse  flow  cases)  'h  =  45°  (pitching  moment  negative  about  Zb 

axis) 


Figure  3.3:  Drag  and  pitching  moment  data  for  regular  and  reverse  flow  cases 


3.5  Surface  flow  visualization 

Surface  flow  visualization  was  performed  using  tufts  attached  to  the  suction  side  of 
the  rotor  blade.  The  tufts  were  2.5  cm  long  pieces  of  yarn  taped  at  2.5  cm  spacing, 
forming  a  grid  over  the  entire  suction  side  of  the  blade.  A  high  definition  camera 
was  used  to  record  the  visualization.  These  experiments  were  performed  at  Re  = 
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Figure  3.4:  Variation  of  chordwise  center  of  pressure  location  with  a  at  various  4/ 

1.7  x  105  and  Re  =  3.0  x  105  at  4/  =  30°,  45°,  —30°,  and  —45°.  The  blade  was  pitched 
through  various  a  in  increments  of  5°.  Reynolds  number  had  little  effect  on  the 
surface  flow  patterns.  However,  large  differences  in  surface  flow  pattern  were  observed 
between  forward  and  backward  yaw  cases.  Hence  only  one  Reynolds  number  case  Re 
=  1.7  x  105  is  discussed  in  detail. 


3.5.1  Forward  yaw  Re  =  1.7  x  10" 


A.  T  =  45° 

Pictures  of  tuft  flow  visualization  taken  at  4/  =  45 

At  a  =  185 


in  Figures  3.6(a)-  3.6(d) 


and  Re  =  1.7  x  105  are  shown 
the  tufts  are  mostly  attached,  but  there 
are  a  few  tufts  at  the  effective  leading  edge  moving  in  a  clockwise,  conical  rotation 
pattern  when  looking  from  upstream.  These  rotating  tufts  are  indicative  of  a 
leading  edge  vortex  beginning  to  form.  There  are  also  some  rotating  tufts  at  the 
blade  tip,  indicating  the  existence  of  a  tip  vortex  interacting  with  the  freestream. 
As  the  blade  is  pitched  to  a  =  190°,  the  leading  edge  vortex  grows  along  the 
span.  The  rapid  growth  of  the  vortex  as  a  increases  from  5°  —  10°  correlates 
well  with  the  large  momentum  deficit  increase  seen  by  Lind  |7i|  for  a  NACA  0012 
airfoil  at  a  =  189°.  In  addition,  there  is  a  small  region  of  the  tufts  exhibiting 
chaotic  motion  near  the  blade  root,  a  characteristic  of  separated  flow,  indicating  a 
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Figure  3.5:  Tuft  visualization  at  Re  =  1.7  x  105 
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vortex  breakdown  region.  The  vortex  breakdown  may  be  attributed  to  the  adverse 
pressure  gradient  caused  by  the  mount  at  the  blade  root.  As  the  blade  is  pitched 
to  a  =  195°,  the  vortex  breakdown  region  grows  larger,  which  correlates  well 
with  the  divergence  of  the  experimental  lift  curve  from  vortex  lift  theory  seen  in 
Figure  3.2(a)  By  a  =  215°,  the  leading  edge  vortex  is  almost  completely  broken 


down  and  there  is  only  a  small  region  of  attached  flow  between  the  remnants  of 
the  leading  edge  vortex  and  the  tip  vortex  region. 

In  the  aerodynamic  loads  section  we  discussed  the  center  of  pressure  on  the  yawed 
blade  in  reverse  flow.  The  physical  reasoning  for  the  variation  of  the  center 
of  pressure  with  angle  of  attack  is  now  reconciled  with  the  flow  visualization 
experiments.  The  growth  of  the  leading  edge  vortex  could  explain  the  center  of 
pressure  movement  seen  in  Figure  3.4(a)  As  the  vortex  becomes  significant  at 
a  =  185°,  the  center  of  pressure  would  shift  toward  the  location  of  the  vortex  at 
the  effective  leading  edge.  As  the  vortex  grows  along  the  span  and  reaches  the 
effective  trailing  edge,  the  center  of  pressure  would  “follow”  the  mean  position 
of  the  leading  edge  vortex  toward  the  effective  trailing  edge,  eventually  settling 
near  the  0.75c  point.  The  vortex  breakdown  region  is  mainly  confined  to  an  area 
near  the  effective  leading  edge  of  the  blade.  The  location  of  the  vortex  breakdown 
region  near  the  effective  leading  edge  could  explain  why  the  chordwise  center  of 
pressure  appeared  to  be  farther  from  the  effective  leading  edge  for  higher  4/  cases 
in  Figure  |3.4(a)  and  Figure  |3.4(b)  The  breakdown  region  would  be  expected  to 


produce  minimal  or  no  lift,  which  would  shift  the  chordwise  center  of  pressure 
toward  the  effective  trailing  edge.  The  lower  T  test  cases  would  not  develop  a 
large  vortex  breakdown  region  because  there  would  be  less  spanwise  flow  and 
a  smaller  adverse  pressure  gradient  at  the  blade  root,  so  the  center  of  pressure 
would  not  be  shifted  as  far  toward  the  effective  trailing  edge. 

It  is  worth  noting  that  the  position  and  shape  of  the  attached  span-wise  vortex 
resembles  that  of  the  lower  surface  vortex  predicted  by  Potsdam  et  al  20  ,  which 


was  attributed  to  ‘reverse-chord  dynamic  stall’.  The  similar  vortex  observed  in 
this  experiment  on  a  static  yawed  blade  suggests  that  the  vortex  predicted  by 
Potsdam  et  al.  may  not  be  an  artifact  of  rotation,  or  pitch  rate.  Based  on  prior 
flow  visualization  experiments  Raghav  et  al.  |:3|]  proposed  a  flow  topology  for  a 
blade  in  forward  yaw  in  reverse  flow  regime. 


B.  T  =  30° 

Pictures  of  tuft  flow  visualization  taken  at  4/  =  30 
in  Figure  3.5(e)  and  Figure  [3.5(f) 


and  Re  =  1.7  x  105  are  shown 
At  a  =  195°  the  observed  flow  pattern  is 
very  similar  to  the  one  seen  at  T  =  45°,  but  at  a  =  215°  a  much  larger  vortex 
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breakdown  region  is  observed.  The  attached  flow  region  is  also  much  smaller  than 
the  one  observed  for  4/  =  45°.  Reasons  for  these  differences  will  be  discussed  in 
later  sections. 


3.5.2  Backward  yaw  Re  =  1.7  x  10" 


A. 


T  =  -45° 

Pictures  of  tuft  flow  visualization  taken  at  4/  =  —45 

At  a  =  195 


in  Figure  3.5(e)  and  Figure  3.5(f) 


B. 


and  Re  =  1.7  x  105  are  shown 
the  tuft  behavior  suggests  that  a 
large  leading  edge  vortex  is  present,  and  no  significant  signs  of  vortex  breakdown 
are  observed.  Unlike  the  forward  yaw  T  =  45°  case,  mount  obstruction  does  not 
impede  the  growth  of  the  leading  edge  vortex  by  creating  an  adverse  pressure 
gradient.  A  smaller  adverse  pressure  gradient  may  be  present  at  the  rotor  blade 
tip  due  to  a  tip  vortex,  but  at  this  angle  of  attack  it  appears  that  the  spanwise 
flow  mostly  negates  the  effect  of  any  tip  vortex.  The  attached  flow  region  is  fairly 
small,  which  is  partly  due  to  mount  interference  near  the  blade  root.  At  a  =  210° 
the  leading  edge  vortex  more  closely  resembles  the  one  seen  at  T  =  45°. 

T  =  -30° 

Pictures  of  tuft  flow  visualization  taken  at  T  =  —30°  and  Re  =  1.7  x  105  are 
shown  in  Figure  3.5(e)  and  Figure  3.5(f)[  At  a  =  195°  a  leading  edge  vortex  is 
apparent,  but  it  does  not  appear  to  have  much  strength.  Some  tufts  near  the 
blade  tip  signify  a  reversal  of  the  expected  span-wise  flow  direction.  At  a  =  215° 
the  rotor  blade  planform  can  be  evenly  divided  into  two  regions:  a  region  where 
the  flow  appears  to  be  fairly  laminar  but  directed  toward  the  root  of  the  blade  in 
opposition  to  the  freestream  flow,  and  a  region  of  chaotic  flow. 


3.5.3  Re  =  3.0  x  105  and  a  =  195° 

The  flow  patterns  for  the  Re  =  3.0  x  105  test  cases  were  essentially  unchanged  from 
the  Re  =  1.7  x  105  cases,  which  is  consistent  with  previous  literature  on  sharp-edged 
delta  wings  |48j. 

3.5.4  Flow  topology 

Raghav  et  al.  |3|  proposed  a  flow  topology  for  a  static  forward  yawed  blade  in  the 
reverse  flow  regime.  The  flow  topology  shows  three  distinct  flow  regimes  -  attached 


CHAPTER  3.  FIXED  WING  RESULTS 


37 


(a)  W  =  45° 


(b)  =  30° 


(c)  4-  =  -45° 


(d)  =  -30° 


Figure  3.6:  Tuft  visualization  at  Re  =  3.0  x  105  and  a  =  195° 


flow,  vortex  flow,  and  separated  flow  demarcated  by  the  dashed  lines  (see  Figure  3.7(a) ). 
In  this  section,  we  present  generalized  flow  topologies  for  a  static  backward  yawed 
blade  in  the  reverse  flow  regime  (see  Figures  3.7(b) -3.7(d) ).  Clearly,  the  surface  flow 
topology  for  a  backward  yawed  blade  differs  significantly  from  that  for  a  forward 
yawed  blade  in  reverse  flow  regime.  Furthermore,  unlike  for  a  forward  yawed  blade 
the  surface  flow  topology  on  a  backward  yawed  blade  also  changes  with  yaw  angle  and 
angle  of  attack.  At  T  =  —30°  the  surface  flow  topology  consists  of  the  familiar  vortex 
and  attached  flow  in  addition  to  regions  of  span-wise  flow  in  direction  opposite  to  the 
expected  span-wise  flow.  In  the  case  of  T  =  —30°  the  flow  topology  also  changes  with 
angle  of  attack  as  observed  from  Figures  3.7(c)  -3.7(d) .  The  effect  of  yaw  angle  of  the 


flow  topology  is  observed  by  comparing  Figure  3.7(b)  to  Figures  3.7(c)  -pk7(d)[  At  the 
time  of  this  writing  about  the  tuft  visualization  it  was  not  clear  what  caused  these 
flow  topology  variations  at  the  various  conditions.  Global  velocity  field  measurements, 
converted  to  static  pressure  distributions,  did  shine  some  light  on  this  behavior  in 
the  rotor  case. 


In  addition,  proposed  flow  topologies  for  a  rotor  blade  in  backward  yaw  in  the  reverse 
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flow  regime  are  depicted  in  Figure |3~.7(b)  -3.7(d) 
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Given  that  the  forward  yaw  flow  topology  regimes  are  fairly  comparable  over  the 
range  of  yaw  angles,  it  is  amenable  to  be  used  for  further  analysis  of  the  physics  of 
the  flow  field.  In  this  section  we  further  analyze  the  variations  in  the  flow  topology 
for  the  forward  yawed  static  blade  at  two  different  yaw  angles  as  the  angle  of  attack 
of  the  blade  changes.  We  first  present  the  variations  in  the  area  occupied  by  the 
different  flow  regimes  with  angle  of  attack  followed  by  the  analysis  of  the  data. 


The  backward  yaw  flow  topology  shows  some  similarities  to  the  forward  yaw  flow 
topology:  there  is  a  region  of  attached  flow  where  the  freestream  flow  first  makes 
contact  with  the  blade  and  a  region  of  vortex  flow  immediately  adjacent  to  it.  A 
key  difference  between  the  two  flow  patterns  is  that  the  backwards  yaw  case  has  a 
region  of  fairly  orderly  flow  bordering  the  far  edge  of  the  vortex  flow  region  while  the 
forward  yaw  topology  features  a  region  of  chaotic,  separated  flow.  In  forward  yaw  the 
blade  mount  forms  a  strong  adverse  pressure  gradient  which  may  accelerate  vortex 
breakdown,  but  in  the  backward  yaw  case  the  leading  edge  vortex  does  not  encounter 
such  an  obstruction  and  is  allowed  to  grow  larger.  It  is  important  to  note  that  the 
flow  visualization  observed  at  4/  =  30°  did  not  follow  the  prescribed  flow  pattern  and 
showed  a  spanwise  flow  direction  that  opposed  the  freestream  flow. 


Figures  3.8(a)  and  |3.8(bJ  illustrate  the  variation  of  the  area  occupied  by  each  flow 
regime  as  a  function  of  angle  of  attack  for  forward  yaw  angles  at  Re  =  1.7  x  105.  An 
insignificant  portion  of  the  blade  area  (<10%)  is  typically  occupied  by  the  tip  vortex 
interaction  region  and  has  not  been  presented  in  the  figures.  The  salient  features  of 
each  flow  regime  are  clearly  observed  in  the  figures:  a )  for  T  =  30°  the  attached  and 
vortex  flow  regimes  monotonously  decrease  with  an  increase  in  angle  of  attack,  while 
the  unattached  flow  increases  with  angle  of  attack,  b)  At  T  =  45°  the  attached  flow 
still  exhibits  a  monotonous  decrease  with  increase  in  angle  of  attack,  however,  the 
vortex  and  unattached  flow  regimes  reach  an  equilibrium  as  angle  of  attack  increases. 


In  order  to  further  understand  the  behavior  of  the  flow  regimes  with  change  in 
yaw  angle  and  angle  of  attack,  we  shall  consider  the  classical  vorticity  transport 
equation  for  an  incompressible,  barot.ropic  fluid  with  conservative  body  forces  [51 1. 
The  following  equations  are  defined  as  per  the  blade  fixed  coordinate  axes  defined  in 


Figure  2.1 


—  =  —  v  ■  Vw  +  u:  ■  Vv  +  vX2w  (3.3) 

where  the  term  representing  the  diffusion  of  vorticity  due  to  viscous  effects,  z/V2u;, 
can  be  neglected  based  on  the  argument  that  the  Reynolds  number  is  sufficiently  high. 
In  addition,  for  simplicity  of  the  argument  if  we  neglect  the  vortex  stretching/tilting 
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Figure  3.8:  Proportion  of  blade  area  occupied  by  the  different  flow  regimes 
function  of  angle  of  attack  at  Re  =  1.7  x  105 


as  a 


term,  uj  ■  V«,  we  arrive  at  a  balance  between  span- wise  vorticity  (ujz 
advection  terms,  such  that: 


production  and 


duj~  du>z  duz  dujz 
+  u— - b  v— - b  w— —  =  (J 


dt 


dx 


dy 


dz 


(3.4) 


From  the  surface  flow  visualization  it  can  be  clearly  observed  that  the  tufts  indicate  a 
predominant  span-wise  velocity  w,  with  negligible  velocity  in  the  Xb  and  yb  directions 
u  and  v  .  The  only  terms  that  remain  in  equation  |3.4| are  the  vorticity  production  term 


and  the  advection  term  in  the  Zb  direction,  such  that  =  0.  Considering 

this  simplified  vorticity  transport  equation,  it  can  be  postulated  that  for  greater  yaw 
angles,  an  equilibrium  state  can  be  achieved  as  long  as  the  vorticity  produced  at  the 
leading-edge  is  rapidly  advected  away  in  the  span- wise  direction.  This  argument  is  now 
used  to  justify  the  equilibrium  state  that  is  achieved  at  the  greater  yaw  angle  4/  =  45° 


as  observed  in  Figure  3.8(b)  when  compared  to  the  lower  yaw  angle  4/  =  30°  as 
observed  in  Figure  3.8(a).  At  4/  =  45°  the  expected  span-wise  velocity  is  greater  (by 
virtue  of  higher  yaw  angle)  than  at  4/  =  30°  and  hence  the  equilibrium  state  is 


easily  achieved.  It  should  be  noted  that  equation  3.3  does  not  include  the  term  that 
accounts  for  destruction  of  vorticity  due  to  annihilation.  The  influence  of  annihilation 
on  the  reverse  flow  aerodynamics  is  yet  to  be  ascertained  from  quantitative  flow  held 
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measurements. 


3.6  Yawed  Rotor  Blade  Held  Stationary 

Velocity  field  measurements  using  stereo  particle  image  velocimetry  are  used  to  exam¬ 
ine  the  sharp-edge  vortex  forming  under  a  rotor  blade  in  reverse  flow.  Comparisons 
with  load  measurements  and  surface  tuft  visualization  show  the  similarities  and 
differences  between  the  rotor  blade  case  and  a  swept  wing  in  reverse  flow.  The  fixed 
wing  case  is  amenable  to  analysis  using  the  Polhamus  suction  analogy,  and  exhibits 
vortex  bursting.  The  rotor  case  shows  the  added  effects  of  radial  acceleration. 

These  are  correlated  with  earlier  fixed-wing  aerodynamic  loads  and  flow  visualization 
studies  on  the  same  rotor  blade  operated  in  forward  and  reverse  facing  yaw  in  the 
wind  tunnel.  The  blade  is  also  held  fixed  at  several  rotor  azimuths  and  angles  of 
attack  to  ascertain  the  rotation  effects  on  the  flowfield.The  circulation  of  the  reverse 
flow  attached  vortex  is  of  the  same  order  of  magnitude  as  the  bound  circulation  of 
the  airfoil  section,  suggesting  that  the  vortex  makes  a  significant  lift  contribution. 

This  section  presents  SPIV  data  for  a  static  yawed  blade  in  reverse  flow  by  exploring 
the  effects  of  reverse  flow  on  rotating  rotor  blades  in  the  John  Harper  2. 13m x  2. 74m 
low  speed  wind  tunnel.  Stereoscopic  particle  image  velocimetry  data  are  used  to 
acquire  three  dimensional  velocity  field  measurements. 


3.6.1  Instrumentation 


All  experiments  were  performed  in  the  2.13mx2.74m  John  Harper  wind  tunnel  at  the 
Georgia  Institute  of  Technology.  The  rotor  blade  used  to  collect  SPIV  data  for  the 
static  experiments  is  an  untwisted,  rectangular  blade  with  a  NACA  0013  planform. 
A  picture  of  the  rotor  blade  used  in  the  static  experiments  is  shown  in  Fig.  |2.1|  See 


Raghav  et  al.  |3|  for  a  detailed  description  of  the  static  yawed  blade  setup.  The 
rotor  used  in  the  rotating  experiments  is  composed  of  two  untwisted,  rectangular, 
NACA  0013  rotor  blades  (identical  to  the  one  used  in  the  static  experiments)  at¬ 
tached  to  a  teetering  rotor  hub  turned  by  a  3.73kW  motor.  A  CAD  model  of  the 


rotor  setup  and  laser  sheet  is  shown  in  Fig.  3.11  and  a  close-up  picture  of  the  rotor 


blade  and  laser  arm  is  shown  in  Fig.  3.10  Rotor  parameters  are  displayed  in  Table  3.3 


Tunnel  speed  was  set  by  gauging  the  dynamic  pressure  and  ambient  pressure  and 
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temperature.  Dynamic  pressure  was  measured  using  a  Barotron  pressure  transducer 
connected  to  a  Pitot-static  probe  in  the  test  section.  The  angular  velocity,  D,  of  the 
rotor  was  measured  using  a  photoelectric  sensor. 

TPP  images  were  collected  using  high-definition  cameras  aligned  with  the  rotor  hub 
plane.  Alignment  was  performed  using  a  laser  balance  and  a  digital  protractor.  A 
camera  was  mounted  on  a  heavy  tripod  two  rotor  diameters  downstream  of  the  rotor 
and  set  at  0°  inclination  using  the  digital  protractor.  The  laser  balance  was  used  to 
set  the  height  and  yaw  of  the  camera  so  that  its  lens  was  directly  pointing  at  the 
center  of  the  rotor  hub.  Another  camera  was  mounted  outside  the  tunnel  test  section 
at  the  rotor  azimuthal  angle  of  4>  =  90°.  This  camera  was  also  aligned  with  the  center 
of  the  rotor  hub  using  the  same  method  as  before.  A  calibration  board  was  used  to 
correlate  pixel  length  with  physical  distance  in  the  plane  of  the  rotor. 

Illumination  for  PIV  was  provided  by  a  Nd:YAG  laser  with  a  pulse  energy  of  200mJ. 
The  laser  was  double  pulsed  to  produce  532nm  visible  light.  A  laser  arm  was  used  to 
deliver  the  beam  to  the  measurement  plane  by  creating  a  light  sheet  with  a  thickness 
of  2mm.  The  flow  was  seeded  with  10pm  droplets  generated  using  a  Laskin-nozzle 
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aerosol  generator.  PIV  images  were  captured  by  a  Imagerlntense  CCD  camera  with 
a  1600  x  1200  pixel  viewing  area.  The  CCD  sensor  had  a  scan  rate  of  16  Hz.  The 
cameras  and  laser  unit  were  synchronized  using  a  programmable  timing  unit. 

Table  3.3:  Rotor  Parameters 


Parameter  Value  Units 


Blade  mass  total 
Blade  span 
Blade  chord 
Blade  aspect  ratio 
Airfoil 
Disc  radius 
Solidity 
Precone 
Motor 
Height 


1.747 

kg 

0.622 

meter 

0.178 

meter 

3.49 

NACA  0013 

0.889 

meter 

0.0895 

1.6 

deg 

3.73 

kW 

1.298 

meter 

3.6.2  Flow  and  Test  Conditions 

The  test  cases  are  summarized  in  Table  13.41  and  described  below: 

1.  The  first  test  collected  images  of  the  TPP  at  ji  =  0.6,  0.8  and  1.0  corresponding 
to  D  =  31. 4r ad/s  and  varying  freestream  speeds.  Images  were  collected  using 
high  definition  cameras  placed  downstream  of  the  rotor  at  T  =  0°  and  to  the 
side  of  the  rotor  at  4/  =  90°.  This  test  established  flow  conditions  that  could  be 
used  for  future  PIV  tests. 

2.  The  second  test  measured  the  chordwise  velocity  field  on  a  static  rotor  blade  in 
reverse  flow  using  PIV.  A  Reynolds  number  of  0.63  x  105  was  chosen  to  match 
the  effective  Reynolds  number  of  a  subset  of  the  rotating  measurements.  A 
test  was  also  performed  at  Re  =  1.70  x  105  for  comparison  to  previous  flow 
visualization  performed  on  a  static  yawed  blade.  Angles  of  attack  were  chosen 
to  match  the  effective  angle  of  attack  of  the  rotating  blade  at  the  same  yaw 
angle. 
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Table  3.4:  Test  Matrix  For  Experiments 


Test 

Data  Acquired 

Conditions 

1 

TPP  Orientation 

li  =  0.6,  0.8, 1.0,  fl  =  31Arad/s 

2 

PIV  Static  Blade 

a  =  191.2°,  192°,  195° 

T  =  30°  Re  =  0.63  x  105, 1.70  x  105 

3 

PIV  Rotating  Blade 

a  =  192°  T  =  240°,  270°, 
li  =  0.8, 1.0  Re  =  0.63  x  105 

3.  The  third  test  acquired  phase  locked  PIV  measurements  on  a  rotating  rotor 
blade  at  advance  ratios  of  //  =  0.8  corresponding  to  Q  =  20.93 rad/s  and 
H  =  1.0  corresponding  to  Q  =  15.7 rad/s.  Measurements  were  made  at  the 
radial  location  of  r/R  =  0.5  and  phase  locked  at  4  =  240°  and  T  =  270° 
equivalent  to  T  =  0°  and  T  =  30°  for  the  static  blade. 


3.6.3  Data  Acquisition/Processing 

TPP  orientation  was  determined  by  measuring  lateral  and  longitudinal  tilt  of  the 
rotor  disc  using  TPP  images.  TPP  tilt  was  measured  by  using  an  image  processor  to 
calculate  the  horizontal  and  vertical  pixel  position  of  the  far  ends  of  the  rotor  hub 
and  performing  an  inverse  tan  operation. 

Davis  7.2  software  by  LaVision  was  used  to  process  the  PIV  data.  Particle  velocities 
were  calculated  from  the  spatial  stereoscopic  cross-correlation  of  the  images.  A  64x64 
pixel  interrogation  window  with  50%  overlap  and  a  second  interrogation  pass  with  a 
reduced  window  size  (32  x  32  pixel)  was  used  to  yield  a  vector  spatial  resolution  of 
1.37mm  which  translates  to  0.0008c.  Post-processing  of  the  vector  images  consisted 
of  applying  a  median  filter  and  light  vector  smoothing  (3  x  3)  for  better  visual 
presentation. 

3.6.4  Averages  flow  field  results 

All  of  the  velocity  vector  fields  presented  are  fixed  in  the  wind  tunnel  coordinate 
system  and  have  not  been  corrected  for  the  blade  pitch  angle.  The  right  hand 
coordinate  system  for  the  vector  field  is  established  so  that  the  positive  y  axis  is  out 
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of  the  plane  and  negative  radial  (spanwise)  velocity  indicates  a  flow  directed  toward 
the  root  of  the  blade.  All  dimensions  are  normalized  by  the  chord  length,  where 
x/c  =  0  corresponds  to  the  geometric  leading  edge  of  the  blade  and  x/c  —  1  is  the 
geometric  trailing  edge  of  the  blade.  The  blue  regions  in  the  velocity  held  figures 
correspond  to  areas  where  PIV  data  were  not  collected. 


Static  Blade 


The  main  motivation  for  performing  PIV  chordwise  velocity  measurements  on  the 
static  blade  in  reverse  how  is  so  that  direct  comparisons  can  be  made  to  measurements 
taken  on  the  rotating  blade.  Therefore,  as  stated  in  the  ‘Flow  andTest  Conditions’ 
section,  T,  a,  and  Reynolds  number  were  chosen  to  match  conditions  present  on  the 
rotating  blade.  One  test  was  also  performed  to  match  conditions  for  the  tuft  how 
visualization  performed  by  Raghav  et  al.  |3|  on  the  same  static  blade  configuration. 
All  static  PIV  tests  were  performed  at  4/  =  30°,  which  corresponds  to  the  blade  being 
yawed  into  the  how  (see  Fig  |2.1|for  the  definition  of  positive  T). 


Figure  3.12  -  3.14  illustrates  the  time  averaged  velocity  vector  overlaid  on  contours  of 
vorticity  at  T  =  30°,  a  =  191.2°  —  195°,  and  Re  =  0.63  x  105  and  1.70  x  I05.  The 
velocity  held  data  taken  at  Re  =  1.70  x  105  and  a  =  195°  corroborate  the  attached 
vortex  seen  using  tuft  visualization  by  Raghav  et  al.  1 3 1  at  the  same  conditions.  All 
of  the  velocity  helds  confirm  the  existence  of  an  attached  vortex  in  reverse  how  over 
a  static  yawed  blade.  Note  that  the  contours  of  vorticity  for  each  velocity  held  have 
different  color  scales,  so  the  strength  of  the  attached  vortex  is  increasing  with  a. 


The  SPIV  data  collected  in  these  experiments  allows  for  interpretation  of  the  radial 
component  of  velocity.  Figure [3T5] shows  the  time  averaged  velocity  vectors  overlaid  on 
the  time  averaged  radial  velocity  component  ( Ur  positive  for  out  of  plane).  The  radial 
component  of  velocity  is  primarily  directed  into  the  plane  of  the  paper,  indicating 
a  spanwise  how  towards  the  blade  root.  This  result  is  expected  for  a  blade  yawed 
into  the  how  and  agrees  with  the  tuft  how  visualization  for  a  static  reverse  how  blade 
presented  by  Raghav  et  al.  |3|.  The  radial  velocity  has  a  lower  magnitude  close  to  the 
surface  of  the  blade  and  increases  further  away  from  the  blade  surface.  This  radial 
velocity  held  can  be  used  to  analyze  the  effect  of  rotation  on  reverse  how,  where  it 
is  expected  to  predominantly  affect  the  radial  velocity  component  due  to  additional 
centrifugal  effects. 
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3.7  Exploratory  Measurements  Under  the  Rotating 
Blade 


In  conjunction  with  the  yawed-blade  tests  above,  a  few  exploratory  tests  were  per¬ 
formed  with  the  rotor  rotating.  While  an  attempt  was  made  to  ’match’  flow  conditions 
as  closely  as  practical  with  the  rotating  blade,  the  comparison  is  very  rough,  as  there 
is  no  spanwise  variation  of  the  relative  freestream  or  the  induced  angle  of  attack. 
Thus  these  experiments  merely  attempted  to  glimpse  what  the  velocimetry  would 
show  in  the  same  stations  were  the  yawed  fixed  blade  was  held. 


Figure  3.16  shows  the  time  averaged  velocity  vectors  overlaid  on  the  time  averaged 
radial  velocity  component  for  a  rotating  rotor  blade  at  /i  =  0.8,  =  20.93rad/s,  and 

r/R  =  0.5  phase  locked  at  T  =  240°  (equivalent  to  a  static  rotor  blade  at  T  =  30°. 
The  radial  velocity  has  a  positive  value  (corresponding  to  flow  directed  form  the 
root  to  the  tip  of  the  blade)  near  the  blade  surface  due  to  reactive  centrifugal  forces 
that  arise  due  to  rotation.  The  radial  velocity  becomes  negative  farther  from  the 
blade  surface,  as  the  radial  stress  driving  the  radial  flow  does  not  extend  far  above 
the  blade  surface.  A  comparison  of  Fig.  3.15  and  Fig.  3.16  demonstrates  the  effect 
of  rotation  on  the  radial  component  of  velocity  in  reverse  flow.  A  comparison  of 
Fig.  3.13  and  Fig.  |3. 17  suggests  that  the  additional  radial  velocity  due  to  rotation 
delays  the  separation  at  the  effective  leading  edge  (geometric  trailing  edge)  at  similar 
flow  conditions. 


With  the  collective  and  cyclic  pitch  settings  of  6 o  =  6°  and  9S  =  —6°,  separation  at 
the  effective  leading  edge  was  not  observed  until  T  =  270°  (equivalent  to  a  static 
blade  at  T  =  0°).  Fig  3.18  shows  an  attached  vortex  on  the  lower  surface  of  the 
blade.  To  clearly  visualize  the  vortex,  the  freestream  velocity  (U0 0)  =  14.88m/s)  was 
subtracted  from  the  velocity  field.  The  circulation  of  the  vortex  (Ty)  was  computed 
using  Stokes  theorem  and  determined  to  be  0.354m2/s.  The  circulation  of  this  vortex 
was  compared  to  the  maximum  bound  circulation  (ra)  of  the  "airfoil"  section  of  the 


wing.  The  bound  circulation  was  computed  using  Eq.  3.5 


r 


a 


1 

2 


£4o  cC Ltjkijz 


(3.5) 


where  c  is  the  chord  length  and  Cmiax  is  the  maximum  static  coefficient  of  lift  taken 
from  aerodynamic  loads  measurement  data  presented  by  Raghav  et  al.  |3|.  The 
computed  bound  circulation  was  0.392m2/s,  which  was  slightly  higher  than  the 
circulation  of  the  vortex  at  the  trailing  edge.  The  fact  that  the  circulation  in  the 
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vortex  is  of  the  same  order  of  magnitude  as  the  bound  circulation  could  imply  that 
the  attached  vortex  which  develops  on  the  lower  surface  of  the  blade  in  reverse  flow 
makes  a  significant  lift  contribution. 

Following  are  the  time  averaged  streamlines  obtained  at  r/R  =  0.5  and  yaw  =  240 
degrees.  The  figures  for  cases  /i  =  0.8  and  1.0  show  rotor  speed  subtracted  stream 
lines  and  trailing  edge  vortex. 

The  last  three  plots  correspond  to  fixed  fi  of  0.85  and  r/R  is  changing  from  0.4  to 
0.5  to  0.6.  The  plots  are  obtained  after  subtracting  the  rotor  cord  wise  velocity  in  X 
direction  and  the  free  stream  velocity  in  Z  (direction  perpendicular  to  the  illuminated 
plane).  Here  we  observe  that  as  we  go  to  higher  r/R,  the  vortex  strength  is  decreasing 
(because  the  reverse  flow  is  decreasing  as  /i  is  constant).  The  other  observation  is 
that  the  in-plane  velocity  in  the  vortex  region  decreases  with  increase  in  r/R  due  to 
centrifugal  effects. 

As  mentioned  before,  these  were  exploratory  measurements.  They  do  show  that  there 
are  vortical  structures  to  be  seen  under  the  rotating  blade. 


Conclusions  from  Fixed  Wing  Experiments 

An  NACA003  rotor  blade  was  first  used  as  a  fixed  wing  in  a  wind  tunnel,  with  air 
loads  measured  at  angle  of  attack,  at  various  forward  and  backward  sweep  angles  in 
both  forward  and  reverse  flow  at  various  Reynolds  numbers.  These  measurements 
suggested  the  presence  of  substantial  amounts  of  vortex  lift.  Lift  curves,  drag  polars, 
and  pitching  moments  were  analyzed  from  low-pass  filtered,  multi-axis  load  cell 
data.  The  Polhamus  suction  analogy  for  delta  wings  was  adapted  for  application 
to  a  yawed  rotor  blade,  and  the  predictions  from  the  mathematical  model  were 
compared  to  experimental  results.  Tuft  visualization  was  performed  at  various  yaw 
angles  and  Reynolds  numbers  to  reconcile  the  aerodynamic  loads  data.  Surface 
tuft  visualization  showed  dramatically  different  vortex  flow  topologies  under  forward 
and  reverse  sweep,  and  forward- facing  versus  reverse-facing  flow.  Vortex  bursting 
appeared  to  be  occurring. 

The  blade  was  then  returned  to  a  2-bladed  teetering  rotor  setup.  It  was  initially  held 
fixed  at  various  azimuths  to  obtain  the  velocity  held  under  the  blade  in  reverse  how 
using  stereo  particle  image  velocimetry.  The  rotor  was  next  operated  at  advance  ratios 
from  0.8  to  1.0.  Velocity  measurements  have  been  completed  at  several  locations  at 
270  degrees  azimuth  and  240  degrees  azimuth.  Other  azimuths  are  to  be  investigated 
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a  well.  The  rotating  rotor  velocity  field  dramatically  confirms  the  formation,  growth 
and  bursting  of  a  sharp-edge  vortex  at  the  240  degrees  azimuth.  The  strong  velocity 
along  the  rotor  radius,  directed  inboard  along  the  core  of  the  vortex  at  240  degrees 
azimuth,  is  reminiscent  of  delta  wing  /  strake  vortices.  It  appears  to  burst,  with  a 
sudden  growth  of  the  core  and  decay  and  spread  of  the  radial  velocity. 

Pitching  moment  data  under  fixed  wing  conditions  confirmed  the  role  of  the  vortex 
in  moving  the  center  or  pressure.  In  the  rotating  flow  case,  the  circulation  in  the 
sharp-edge  vortex  is  again  found  to  be  substantial  and  of  the  same  order  as  the  overall 
bound  circulation  over  the  airfoil  section,  and  thus  bound  to  cause  first-order  effects 
on  the  pitching  moment  of  the  blade.  The  results  have  thus  conclusively  proven  the 
existence  of  strong  vortex  lift  characteristics  in  the  aerodynamics  of  a  blade  in  reverse 
flow.  The  conclusions  are  enumerated  in  summary  below: 

1.  The  reverse  flow  lift  curve  slopes  were  higher  than  the  corresponding  regular 
flow  lift  curve  slopes  at  forward  yaw  angles.  The  reverse  flow  lift  curve  slopes 
at  backward  yaw  angles  were  lower  than  the  reverse  and  regular  flow  lift  curve 
slopes  at  forward  yaw  angles  and  did  not  follow  the  expected  cos2T  fit.  The 
observed  flow  behavior  of  the  backward  yaw  cases  also  differed  significantly  from 
the  forward  yaw  cases.  The  lift  asymmetry  between  the  forward  and  backward 
yaw  angles  could  affect  trim  conditions  in  very  high  advance  ratio  flight.  It  is 
important  to  note  that  mount  interference  could  have  influenced  test  results  for 
the  backward  yaw  angle  cases,  which  represent  a  helicopter  blade  in  the  rotor 
azimuth  regime  T  >  270°. 

2.  The  lift  predictions  made  by  the  Polhamus  delta  wing  equations  showed  excellent 
correlation  with  experimental  results  at  forward  yaw  angles  of  T  >  30°.  The 
Polhamus  equations  showed  fair  correlation  at  0°  <  T  <  30°,  which  is  to  be 
expected.  The  ability  of  the  Polhamus  equation  to  accurately  predict  static 
loads  in  the  reverse  flow  regime  set  the  stage  to  separate  out  the  effects  of 
rotation,  from  experiments  on  a  rotating  blade. 

3.  Drag  polars  showed  that  the  regular  flow  regime  produced  higher  L/D  values 
than  the  reverse  flow  regime  for  the  yawed  rotor  blade.  Pitching  moment  data 
suggests  that  the  chordwise  center  of  pressure  is  near  0.75c  in  the  reverse  flow 
regime.  The  larger  aerodynamic  loads  in  the  reverse  flow  region  (as  evidenced 
by  the  higher  lift  curve  slopes)  coupled  with  this  shifted  center  of  pressure 
would  result  in  high  torsional  loads  and  pitch  link  loads  during  flight  in  the 
reverse  flow  regime. 
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4.  Surface  tuft  flow  visualization  showed  strong  evidence  of  an  attached  leading 
edge  vortex  in  the  reverse  flow  regime.  The  vortex  size  was  strongly  dependent 
on  the  angle  of  attack  and  correlated  well  to  the  aerodynamic  loads  observed  at 
the  same  yaw  angle.  The  attached  vortex  contributes  to  higher  blade  loading  in 
the  reverse  flow  regime.  The  position  and  shape  of  the  attached  vortex  closely 
resembled  the  “reverse-chord  dynamic  stall”  vortex  predicted  by  CFD  for  the 
slowed  UH-60  rotor  test  1 20 1 .  This  suggests  that  the  vortex  may  not  be  an 
artifact  of  rotation. 

5.  The  flow  topology  over  the  forward  yaw  blade  in  reverse  flow  has  three  main 
regimes:  the  vortex,  attached  and  separated  flow.  The  flow  topology  variation 
with  yaw  angle  suggests  that  at  greater  yaw  angle  configuration  the  size  of  the 
attached  span- wise  vortex  attains  a  maximum  equilibrium  value  when  compared 
to  the  lower  yaw  angle  configuration.  This  observation  is  justified  by  using  a 
vorticity  transport  analysis  of  the  flow  field. 

3.7.1  Implications 

The  fixed  wing  experiments  thus  showed  that  this  rotor  blade  would  develop  a 
sharp-edge  vortex  under  reverse  flow  conditions.  The  blade  aspect  ratio  is  low,  and 
this  raises  questions  about  the  validity  of  the  conclusions  for  helicopter  blades  that 
have  much  higher  aspect  ratio.  However,  these  concerns  are  not  justified,  when  one 
considers  the  nature  of  the  rotor  reverse  flow  problem  in  detail.  For  one  thing,  aspect 
ratio  is  not  the  primary  issue  here.  The  same  blade  could  have  been  built  with  a 
much  smaller  aspect  ratio,  simply  by  reducing  chord.  That  would  not  have  changed 
the  conclusion  that  the  sharp-edge  vortex  would  form;  it  would  only  have  reduced 
the  Reynolds  number  based  on  chord.  Secondly,  consider  that  in  the  case  of  a  rotor, 
the  region  of  reverse  flow  is  generally  limited  to  a  part  of  the  blade  span.  The  highest 
reverse  flow  dynamic  pressure  is  at  the  root.  The  point  of  zero  chordwise  relative 
velocity  reaches  the  tip,  at  advance  ratios  less  than  1,  only  at  270  degrees  azimuth.  So 
again,  the  region  of  reverse  flow  vortex  is  limited  in  practice.  Even  in  full-scale  flight, 
the  reverse  flow  Mach  number  is  small,  so  that  incompressible-regime  tests  are  entirely 
appropriate.  What  is  the  effect  of  the  larger  Reynolds  number?  Probably  a  tighter 
vortex,  but  if  the  aspect  ratio  is  large,  this  vortex  might  again  cover  most  of  the 
blade,  as  is  the  case  here.  Were  a  small-chord  blade  to  be  used  in  these  experiments, 
the  larger-core,  lower  Reynolds  number  vortex  might  have  spilled  over  the  chord,  and 
that  would  be  unrealistic.  Thus  we  argue  that  the  choice  of  experimental  parameters 
is  fully  justified  by  the  results. 
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These  experiments  have  thus  proved  the  first  part  of  the  hypothesis:  That  a  blade  in 
reverse  flow  generates  a  sharp-edge  vortex,  much  like  the  sharp  edge  of  a  delta  wing 
at  angle  of  attack.  The  experiments  done  with  the  blade  held  at  expected  rotor  blade 
attitude,  show  that  there  is  significant  circulation  in  the  vortex.  Polhamus’  suction 
analogy  can  be  extended  to  calculate  the  vortex  lift  and  the  shift  in  center  of  pressure 
due  to  this  vortex,  and  they  are  both  significant.  Now  we  can  move  to  the  next  step. 
What  will  happen  when  the  rotor  blade  is  rotating?  What  are  the  special  effects  of 
rotation?  Given  our  prior  experience  with  studying  the  strong  radial  flow  on  rotor 
blade  on  the  retreating  blade  side,  radial  flow  effects  were  certainly  uppermost  in  our 
minds  as  we  ventured  into  this  regime. 
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Figure  3.10:  Photo  of  the  Stereo  PIV  setup 
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Figure  3.14:  Static  Blade  at  d'  =  30°,  a  =  195°,  Re  =  1.70  x  105 
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Figure  3.15:  Radial  flow  on  static  blade  at  4/  =  30°,  a  =  192°,  Re  =  0.63  x  105 


x/c 


Figure  3.16:  Radial  flow  under  rotating  blade  at  fi  =  Os,  yaw=240,  r/R=0.5 
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Figure  3.17:  Vorticity  Contours  Under  Rotating  Blade  At  ^  =  240°,  /i  =  0.8,  and 
=  20.93rad/s  which  yields  a  =  192°  and  Re  =  0.63  x  105  at  r/R  =  0.5 


x/c 


Figure  3.18:  Vorticity  Contours  Under  Rotating  Blade  At  'h  =  270°,  /i  =  0.8,  and 
hi  =  20.93rad/s  which  yieldscc  =  196°  and  Re  =  0.635  x  105  at  r/R  =  0.5 
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Figure  3.20:  Yaw  =  240,  r/R  =  0.5,  //  =  0.8,  Local  Reynolds  No. 
Observed  approximate  circulation  at  vortex  =  0.75m2/s 
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Figure  3.21:  Yaw  =  240,  r/R  =  0.5,  //  =  1.0,  Local  Reynolds  No.  =  1.89  *  105, 
Observed  approximate  circulation  at  vortex  =  1.44m2/s 
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(a)  Yaw  =  240,  r/R  =  0.4,  fj,  =  0.85,  Rotor  (b)  Yaw  =  240,  r/R  =  0.5,  =  0.85,  Rotor  Speed 

Speed  =  200  RPM  =  200  RPM 


X  (mm) 


(c)  Yaw  =  240,  r/R  =  0.6,  n  =  0.85,  Rotor  Speed 
=  200  RPM 


Figure  3.22:  Radial  variation  of  velocity  field  at  a  fixed  advance  ratio  of  0.85,  at  240 
degrees  azimuth 


(m/s) 


Chapter  4 

Rotor  Blade  Results 


4.1  Averaged  flow  field  data 

SPIV  data  presented  here  are  averaged  over  100  instances.  All  the  streamline  plots 
are  fixed  in  the  blade  reference  frame  by  subtracting  the  blade  velocity  from  the 
averaged  data.  The  contour  plot  in  the  background  of  the  streamline  figures  depicts 
the  out-of-plane  velocity  component  (Ur).  For  clarity,  the  contribution  of  the  free 
stream  (due  to  yaw)  in  the  out-of-plane  direction  has  been  removed  from  the  Ur 
contour  plot.  The  following  equations  show  how  the  velocities  in  the  plots  have  been 
modified. 


Ux, Blade  ~  Ux,Lab  ~  fxQ 


(4.1) 


Ur, contour  Ur^0i)ServecL 


U freestream%£> 'i'H' (^70 


(4.2) 


where  U freestream  is  the  wind  tunnel  flow  speed.  The  right  hand  coordinate  system 
has  been  used  in  all  plots,  so  the  positive  r  direction  points  out  of  the  plot.  Positive  X 
direction  points  towards  leading  edge  of  blade.  The  region  of  streamline  plots  masked 
in  blue  represents  the  blade  geometry  and  that  masked  in  gray  represents  the  region 
where  data  were  not  obtained  due  to  laser  shadow.  The  laser  sheet,  and  thus  the 
images,  cut  the  blade  perpendicularly  in  all  cases. 
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4.1.1  Azimuth  240  deg. 


At  0  =  240,  the  blade  is  in  a  30-degree  forward  yaw  condition.  Figure  44]  displays 
streamline  plots  at  different  sections  along  the  blade  at  the  fi  =  0.85  condition.  All 
the  plots  (except  that  at  r/R  =  0.7)  have  a  high  strength  vortex  attached  to  the 
blade  trailing  edge.  This  closely  resembles  a  forward  swept,  sharp  edged  wing  at 
angle  of  attack.  The  vortex  strength  and  into-the-plane  velocity  (~Ur)  at  the  core  of 
the  vortex  decrease  with  increasing  r/R  (moving  towards  the  tip).  One  reason  for  this 
behavior  is  the  decrease  in  reverse  flow  going  from  hub  to  tip  due  to  the  increase  in 
tangential  speed  with  r/R.  The  streamline  plot  at  0.7  r/R  does  not  show  any  vortex 
because  of  negligible  reverse  flow  at  that  section. 


Figure  4.1:  Streamline  plots  from  SPIV  near  the  sharp  edge  of  the  blade  at  different 
r/R  at  0  =240  and  /i  =0.85. 


The  increase  in  vortex  strength  parallels  the  rise  in  reverse  flow  with  increase  in  /i 
as  apparent  from  Figure  4.2  There  is  almost  no  reverse  flow  at  /j,  of  0.7  at  r/R  = 
0.6.  High  reverse  flow  near  the  root  causes  a  tight,  strong  vortex.  There  is  no  vortex 
shedding  observed  at  this  azimuth. 


At  azimuth  of  270  deg,  the  blade  is  at  zero  yaw.  In  Figure  |4~j]  All  the  sections  feature 
a  large  vortical  structure  which  is  not  centered  near  the  trailing  edge,  unlike  the  cases 
at  0  =  240  deg.  This  structure  is  considerably  clearer  than  what  one  might  expect 
due  to  simple  recirculation  in  separated  flow.  We  can  postulate  that  separated  flow 
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Figure  4.2:  Streamline  plots  at  r/R  of  0.6  at  0  =240  and  /j  =0.85. 


does  contribute  to  the  strength  of  the  vortex  as  it  hits  the  vortex  tangentially.  The 
vortex  shape  is  more  well-defined  at  sections  near  the  tip  while  they  are  not  so  near 
the  hub.  This  vortex  shape  deformation  is  possibly  due  to  high  reverse  flow  strength 
in  the  region.  Reverse  flow  in  this  case  does  not  cause  significant  increase  in  vortex 
strength  like  in  the  case  of  -0  =  240  deg.  as  the  blade  is  not  in  a  yawed  condition. 


In  Figure  4.4  with  increase  in  /i,  the  vortex  is  further  away  from  trailing  edge  in  the 
flow  direction.  The  streamline  plots  at  /i  of  0.7  and  0.85  with  the  vortex  extending 
well  below  the  blade  edge,  clearly  refute  models  based  on  shedding  from  the  blade. 
At  /i  =  1,  the  vortical  structure  at  r/R  =  0.6  appears  to  be  strong  but  confined  very 
close  to  the  blade  surface.  The  overall  finding  of  the  vortex  well  below  the  blade  edge 
suggests  that  the  vortex  is  created  at  some  0  less  than  270  deg.  where  the  blade  was 
in  forward  yaw  position.  The  vortex  is  then  convecting  with  the  blade.  The  other 
observation  is  that  the  maximum  negative  radial  velocity  (Ur)  is  found  to  be  at  the 
location  where  the  vortex  touches  the  bottom  of  the  blade  for  all  cases  at  this  0.  This 
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Figure  4.3:  Streamline  plots  from  SPIV  near  the  sharp  edge  of  the  blade  at  different 
r/R  at  0  =270  and  n  =0.85. 


can  be  explained  by  considering  the  radial  static  pressure  gradient  along  the  radial 
direction:  the  vortex  being  stronger  inboard,  the  static  pressure  near  the  surface  is 
expected  to  be  lower,  so  that  there  is  a  favorable  gradient  directed  radially  inboard. 
This  would  induce  a  radial  velocity  directed  inboard,  at  the  outer  r/R  locations,  as 
seen.  This  postulate  requires  confirmation  from  static  pressure  data. 


4.1.2  Azimuth  300  degrees 

At  azimuth  of  300  deg.,  the  blade  is  in  a  30  degree  backward-  yaw  condition.  The 
streamline  plots  in  Figure  |4R] resemble  those  at  0  =  270  deg.  in  terms  of  vortex  shape, 
size  and  trends.  Here  too  the  vortex  near  the  hub  appears  less  distinct.  Increasing  // 
leads  to  a  trend  similar  to  that  at  0  =  270  deg.  and  the  vortex  loses  its  shape  with 
increase  in  reverse  flow  (Figure  [4. 4[).  But  this  time,  this  effect  is  seen  to  occur  for 
lower  /i  values.  This  indicates  that  the  large  vortical  structures  observed  at  0  =  270 
deg.  are  still  present  when  the  blade  comes  to  0  =  300  deg.  However,  they  are  not 
as  strong  now  and  are  diffusing  energy.  In  this  case  too,  maximum  negative  values  or 
Ur  are  found  where  the  vortex  interacts  with  the  blade  bottom  surface. 
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4.2  Pressure  Coefficient  Contours 


The  velocity  field  measurements  obtained  from  PIV  data  were  used  to  find  pressure 
coefficients  in  the  flow  field.  Reference  1 52 1  gives  the  details  of  the  process  used. 
Figure  4.7  shows  a  pressure  coefficient  contour  plot  for  the  case  of  i/j  =  2700,  r/R 
=  0.4  and  /i  =  0.85.  Because  the  flow  here  is  primarily  incompressible  and  viscous 
effects  are  restricted  to  boundary  and  vortex  core  only,  the  Cp  found  using  Bernoulli 
equation  is  not  only  an  upper  bound  but  also  a  fairly  accurate  estimate  for  most  of 
the  region. 


The  high  values  of  negative  Cp  seen,  indicate  that  the  vortex  is  highly  three- 
dimensional,  having  grown  along  its  core,  and  must  have  been  generated  at  an 
azimuth  well  ahead  of  270  degrees.  The  highest  negative  values  are  seen  near  the 
surface,  where  the  edge  of  the  vortex  core  touches  the  blade.  Figure  |4~8]  is  the  contour 
plot  of  Cp  at  the  blade  surface  along  the  blade  at  t/j  =  270  deg.  obtained  by  inter¬ 
polating  the  Cp  distribution  along  blade  surface  at  each  section.  The  interpolation 
is  highly  jagged  because  of  the  spare  number  of  planes  where  SPIV  data  have  been 
obtained.  It  shows  that  the  highest  negative  surface  pressures  occur  around  x/c  of 
0.15  from  the  TE  (x/c  of  0.85  measured  from  the  rounded  edge  of  the  blade),  with 
this  region  moving  diagonally  out  to  the  sharp  edge  at  outboard  locations.  Again, 
this  is  consistent  with  the  idea  of  the  vortex  having  been  formed  along  the  sharp 
edge  at  an  earlier  azimuth,  but  convecting  faster  than  the  blade  at  the  more  inboard 
locations. 


Figure  T9  follows  the  discussion  on  pressure  loss  measured  by  Couette  flow  analogy 
as  explained  in  1 52 1 .  It  is  observed  that  pressure  loss  as  compared  to  freestream 
dynamic  pressure  is  of  the  order  of  10~4.  The  seemingly  low  values  of  pressure  losses 
validate  the  use  of  Bernoulli  equation  for  initial  best  guess  and  accurate  estimate 
for  pressure  field.  Pressure  losses  along  streamline  near  to  vortex,  measured  as  an 
artifact  of  shear  stress,  attain  negative  values  in  the  vortex  influenced  regions  and 
goes  to  zero  as  the  streamline  leaves  the  vortex.  As  a  next  step,  the  pressure  loss 
may  be  refined  further  and  can  be  used  as  a  correction  term  to  accurately  predict  the 
static  pressures  in  the  entire  flow  field. 


4.3  Interpretation  of  Results 

Figure  1.2  provides  a  simplistic  illustration  of  the  different  predictive  models  seen  in 
the  literature.  Due  to  probably  oversimplification,  we  are  not  associating  these  with 
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specific  references.  The  most  basic  method  was  to  take  airfoil  data  from  2D  wind 
tunnel  tests  in  the  region  around  180  degrees  angle  of  attack,  and  correct  the  data 
for  yaw  with  a  sinusoidal  correction,  and  if  needed  an  aspect  ratio  correction.  At  the 
next  level,  with  experiments  revealing  evidence  of  vortices  in  the  flowfield,  vortex 
shedding  was  postulated;  both  from  the  blunt  edge  when  using  NACA0012  type 
airfoils  with  blunt  leading  edges,  and  from  the  sharp  edge  as  from  a  thin  flat  plate 
held  at  negative  angle  of  attack.  Karman  vortex  street  analogues  were  also  developed. 
These  models  have  recently  been  adapted  to  include  unsteady  pitching  effects,  again 
from  airfoil  experimental  data.  Other  studies  with  unswept  airfoils  have  modeled  the 
sharp-edge  flow  as  a  diffused,  separated  shear  layer,  with  the  blunt  edge  shedding 
vortices.  Others  have  postulated  a  reverse  dynamic  stall  process,  explaining  the 
occurrence  of  sharp  pitching  moments  as  well  as  strong  vortices  seen  in  the  flowfield 
in  experiments.  At  the  bottom  is  our  postulate:  that  the  flow  around  and  under  a 
sharp,  yawed  edge  behaves  similarly  to  that  around  and  over  a  sharp-edge  delta  wing 
or  other  highly  swept  wing  at  angle  of  attack.  A  strong,  coherent,  helical  sharp-edged 
vortex  is  formed.  Where  this  encounters  an  adverse  pressure  gradient,  as  near  a  hub, 
it  can  burst  and  lift  away  from  the  surface.  The  vortex-induced  lift  and  pitching 
moment  can  be  computed  using  Polhamus  suction  analogy,  as  well  as  models  for 
vortex  flows  at  higher  angles  of  attack. 


Comparing  images  from  azimuths  240  (30  degrees  forward  yaw)  with  those  from 
azimuth  300  (30  degrees  backward  yaw)  reveals  some  of  the  effects  of  rotation  and 
history.  Figure  4. 10|  illustrates  the  difference  in  the  onset  velocity  field  between  the 
fixed,  yawed  wing  and  the  rotating  wing.  The  vortex  strength  increases  as  one  moves 
along  the  blade  span.  The  first  cause  is  the  variation  of  freestream  velocity  and  hence 
dynamic  pressure.  This  effect  implies  a  higher  strength  near  the  root  in  both  cases, 
since  the  reverse  flow  dynamic  pressure  is  highest  at  the  root.  In  analysis,  this  effect 
can  be  easily  included  to  a  model  based  on  2-D  airfoil  aerodynamics. 


The  second  effect  is  the  OintegratingO  effect  along  the  blade  span.  This  is  a  nonlinear 
effect  of  yaw,  and  is  related  to  the  flow  around  the  leading  edge  of  a  delta  wing  or 
high  swept  wing.  At  any  given  section  of  such  a  wing,  the  vortex  observed  is  due  to 
vorticity  generated  by  that  section  as  well  as  the  vorticity  in  the  shear  layers  rolling 
up  over  the  edge  at  sections  upstream  of  the  one  being  observed.  Thus  the  attached 
vortex  observed  at  a  section  can  be  thought  of  as  due  to  OintegratedO  contribution 
by  all  the  upstream  sections.  The  vortex  strength  should  this  increase  along  the  swept 
(yawed)  edge.  This  effect  implies  strength  growing  inboard  in  the  case  of  the  240 
degree  (forward  swept)  case  and  towards  the  tip  in  the  300  degree  (backward-swept) 


case. 
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In  |2|  we  showed  the  following  points  about  the  nature  of  the  flow  field  around  a  rotor 
blade  in  the  reverse  flow  region. 

1.  The  reverse-flow  rotor  blade  encounters  a  significantly  different  onset  velocity 
field  compared  to  a  fixed,  yawed  wing. 

2.  There  is  a  systematic  difference  in  dynamic  pressure  of  the  inset  flow  from  root 
to  tip,  the  root  encountering  the  highest  reverse  flow  speed. 

3.  At  240  degrees  azimuth,  a  strong  vortex  is  generated  along  the  sharp  edge,  and 
grows  in  the  inboard  direction. 

4.  The  increase  in  vortex  strength  directed  inboard  is  due  both  to  the  yaw  effect, 
where  a  helical  vortex  grows  along  the  edge,  and  due  to  the  linear  increase  in 
onset  velocity  moving  inboard. 

5.  The  radial  velocity  is  directed  towards  the  hub  due  to  the  helical  vortex  and 
due  to  the  favorable  pressure  gradient  at  the  surface  along  the  inward  radial 
direction. 

6.  At  270  degrees  azimuth,  the  vortices  seen  are  those  generated  at  azimuths  ahead 
of  270.  The  vortex  is  detached  from  the  blade  and  is  converting  along  with  the 
blade.  The  core  shows  the  effects  of  a  burst  vortex. 

7.  At  300  degrees  azimuth,  the  vortex  structure  looks  similar  to  that  at  270  degrees. 

8.  From  the  above,  a  consistent  model  of  the  blade  aerodynamics  may  be  developed, 
based  on  vortex-flow  concepts  adapted  from  delta  wings  at  angle  of  attack. 

9.  Models  based  on  2-D  airfoil  aerodynamics  and  steady  of  unsteady  shedding,  do 
not  usefully  capture  the  physics  of  the  reverse-flow  rotor  problem. 

10.  Static  pressure  contours  in  the  flow  field  show  the  low-pressure  region  developed 
close  to  the  surface  by  vortex  interaction,  resulting  in  vortex  lift. 

11.  Static  pressure  contours  on  the  blade  surface  at  270  degrees  and  beyond,  show 
the  vortex  convected  onto  the  blade,  with  minimum  pressure  occurring  0.15c 
from  the  sharp  edge  (0.85c  from  the  blunt  edge).  This  supports  the  model  of 
the  vortex  detaching  from  the  edge  as  yaw  becomes  small,  bursting,  and  having 
the  core  edge  flow  come  close  to  the  blade  surface. 

In  continuing  work,  we  explore  the  static  pressure  field. 
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1.  The  pressure  field  computed  from  the  velocity  field  assuming  constant  stagnation 
pressure,  shows  why  there  is  an  inboard-directed  radial  flow  at  270  and  300 
degrees  azimuth. 

2.  The  effect  of  shear  stress  on  the  pressure  loss  is  estimated  along  a  few  streamlines. 
This  shows  that  the  static  pressure  losses  are  dominant  in  vortex  influence  region 

3.  Accordingly  the  pressure  gradient  along  the  blade  surface  is  slightly  lower  than 
the  constant  stagnation  pressure  case. 

4.  Pressure  losses  that  are  influenced  by  shear  stress  near  to  vortex  is  significantly 
low,  thus  validating  the  application  of  Bernoulli  equation  for  nearly  constant 
stagnation  pressure. 


4.4  Surface  Pressure  Contours 


SPIV  data  at  5  radial  locations  starting  from  root  r/R  0.4  to  near  tip  r/R  0.7  were 
subjected  to  an  interpolation  to  satisfy  the  differential  mass  conservation  equation, 
and  the  results  were  interpolated  using  a  cubic  spline.  The  correspondence  with  the 
continuity  equation  is  worst  near  the  vortex  core.  Figure  |2^b)  shows  the  orientation. 
Color  bars  in  Figure  4.11  denote  the  pressure  coefficient  Cp  normalized  by  freestream 
dynamic  pressure.  The  term  "LE"  refers  to  the  blunt  edge  of  the  blade.  Near  the 
tip  (r/R  above  0.85),  there  is  flow  acceleration  ahead  of  the  sharp  edge,  with  some 
suction  for  the  first  10%  of  chord  from  the  sharp  edge  (the  relative  velocity  is  lowest 
at  the  tip).  Strongly  negative  Cp  at  the  root  is  explained  by  the  higher  relative 
reverse  speed.  This  pressure  gradient  negates  centrifugal  effects.  As  /x  increases  from 
0.7  to  1.0,  this  root  to  tip  pressure  gradient  rises. 


Figure  |4T2j  shows  the  situation  at  270  degrees  azimuth.  The  blade  is  perpendicular  to 
the  flow,  but  moving  rapidly  and  with  a  linear  variation  of  onset  velocity  from  root  to 
tip.  History  effects  must  be  considered.  The  very  strong  negative  Cp  regions  occurring 
upstream  of  the  sharp  edge  clearly  indicate  that  a  strong  vortex  is  convecting  from 
upstream  of  the  blade  rather  than  being  caused  by  any  ’shedding’  from  the  blade 
as  2D  models  might  predict.  The  negative  Cp  near  TE  at  r/R  0.4-0.45  raises  the 
possibility  of  flow  separation  or  a  secondary  vortex  being  formed.  The  negative  Cp 
observed  further  downstream  is  certainly  due  to  a  vortex  convecting  with  the  rotor. 
At  azimuth  300,  (Figure  [4. 13[)  the  blade  in  reverse  flow  behaves  like  a  backward  swept 
wing.  Pressure  contours  show  evidence  of  the  vortex  at  lower  r/R.  At  fi  of  1.0  the 
vortex  appears  to  have  convected  further  downstream  and  left  the  image  frame.  The 
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upstream  flow  is  inclined  towards  the  tip  of  rotor,  while  in  the  reverse  flow  region, 
flow  moves  towards  the  root  based  on  the  pressure  gradient.  Being  generated  from  the 
velocity  fields  obtained  at  a  few  discrete  radial  stations,  these  plots  are  still  jagged. 


4.5  Vortex  core  profile 


To  get  beyond  slender  wing  theory  at  small  angle  of  attack,  the  effect  of  flow  detach¬ 
ment  into  first  an  attached  vortex  and  then  a  detached  vortex,  must  be  considered. 
The  velocity  profile  across  the  vortex  was  measured  as  shown  in  Figure  [4.14  at  270 
degrees  azimuth,  a  clearly  ’nonlinear’  case  in  aerodynamics  where  the  vortex  core 
center  is  significantly  away  from  the  surface.  The  result  shows  the  signature  of  a 
burst  vortex:  the  solid-body  core  profile  extends  over  the  full  diameter  of  the  vortical 
flow. 


4.6  Blunt  Edge  Flowfield 


At  the  final  stage  of  the  project,  we  completed  SPIV  data  acquisition  around  the  blunt 
edge  of  the  blade  at  240  and  270  degrees.  The  results  are  qualitatively  in  Figure  4.15 
and  presented  below  in  Figure  |4.16|  The  phase-resolved,  ensemble-averaged  flowfield 


at  r/R  =  0.5  at  240  degrees  azimuth  is  not  massively  separated.  It  in  fact  shows  the 
streamline  deflecting  upwards  around  the  blunt  edge,  resulting  in  a  net  upward  flow 
turning  angle.  This  is  in  accordance  with  expectations,  showing  a  negative  lift  at 
this  azimuth  and  r/R  location.  There  was  considerable  cycle-to-cycle  difference  in 
the  blunt  edge  flowfield,  and  some  images  did  show  separation  or  downward-tilting 
streamlines.  Thus  a  large  fluctuating  lift  and  drag  load  is  expected  due  to  the  nature  of 
the  flow  around  the  blunt  edge.  At  270  degrees  net  streamline  curvature  is  downward, 
suggesting  positive  lift,  generation.  Again  there  was  some  cycle-to-cycle  uncertainty. 
In  all  cases  there  was  strong  but  fluctuating  radial  flow  in  the  very  near  wake,  showing 
that  the  flowfield  around  the  blunt  edge  is  highly  3-D  and  vortical.  Some  of  this  is  to 
be  expected.  The  radial  flow  is  outward  everywhere  other  than  in  the  SEV  zone. 
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4.7  Aspect  Ratio  Effects 

The  aspect  ratio  of  the  blade  used  here  is  small  by  helicopter  rotor  standards.  It 
is  only  3.49  for  the  hard-surfaced  blade  itself,  though  with  root  cut-out  and  hub 
diameter  and  the  opposite  side  included,  the  aspect  ratio  in  airplane  terms  is  well 
above  8.  Still,  values  for  flying  helicopters  are  on  the  order  of  40  or  more.  How 
well  do  the  results  here  capture  expected  phenomena  in  full-scale  flight?  A  thought 
experiment  is  in  order.  The  wind  tunnel  used  here  is  large  by  most  standards:  a  2.74m 
test  section  width  being  the  relevant  dimension.  The  rotor  diameter  is  constrained 
by  wall  effect  avoidance.  The  chord  Reynolds  number  of  the  experiment  is  low,  but 
not  hugely  lower  than  flight  conditions  because  in  reverse  flow,  at  least  one  section  of 
the  blade  goes  through  zero  relative  velocity.  So  the  vortex  structures  seen  here  are 
not  very  unrepresentative  of  flight  conditions.  The  full-flight  vortex  is  expected  to  be 
significantly  tighter  and  stronger  because  of  higher  flow  velocities,  and  hence  perhaps 
occupy  less  of  the  chord  than  is  seen  in  our  experiments. 

The  thought  experiment  is  to  increase  the  aspect  ratio  of  the  present  experimental 
blade  by  simply  using  a  smaller-chord  blade.  What  might  happen?  The  vortex 
generation  process  is  not  expected  to  change  significantly,  since  it  occurs  at  the  sharp 
edge.  The  fact  that  the  blunt  edge  is  encountered  quicker,  will  be  mitigated  by  the 
tighter  vortex  if  we  run  at  flight  conditions.  So  in  net  effect,  the  results  seen  will  not 
be  significantly  different.  Multiple  bursts  and  regeneration  of  the  SEV  may  occur, 
but  this  is  not  seen  on  large  delta  wings  which  have  very  long  leading  edges.  So,  in 
net  effect,  the  vortex  should  not  be  very  different.  Nor  should  the  conclusions  about 
surface  pressure  gradients.  Nor  about  the  3-D  flow  around  the  blunt  edge. 
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Figure  4.4:  Streamline  plots  at  r/R  of  0.6  at  ^  =270  and  /j  =0.85. 
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Figure  4.5:  Streamline  plots  from  SPIV  near  the  sharp  edge  of  the  blade  at  different 
r/R  at  -0  =300  and  /i  =0.85. 


CHAPTER  4.  ROTOR  BLADE  RESULTS 


72 


-7.5  -5.0  -2.5  0.0  2.5  5.0 

Ur  [m/s] 


Figure  4.6:  Streamline  plots  at  r/R  of  0.6  at  ^  —  300  and  /j  =0.85. 
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Figure  4.7:  Contours  of  static  pressure  coefficient  in  the  flow  field  generated  from  the 
SPIV  data  by  assuming  constant  stagnation  pressure.  Azimuth  240,  r/R=0.5,  ^=0.85 
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Figure  4.8:  Contours  of  static  pressure  coefficient  on  the  blade  generated  from  the 
SPIV  data  by  assuming  constant  stagnation  pressure.  Azimuth  270,  ^=0.85 
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Figure  4.9:  Viscous  shear  stress  coefficient  along  a  streamline  through  the  sharp-edge 
vortex,  as  a  measure  of  the  loss  in  stagnation  pressure. 
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Figure  4.10:  Various  models  for  reversed  flow  on  rotor  blades. 
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Figure  4.11:  Static  pressure  contours  obtained  from  the  interpolated  velocity  field,  at 
-0  =  240  for  different  advance  ratios 
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Figure  4.12:  Static  pressure  contours  obtained  from  the  interpolated  velocity  field,  at 
-0  =  270  for  different  advance  ratios 
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Figure  4.13:  Static  pressure  contours  obtained  from  the  interpolated  velocity  field,  at 
if}  =  300  for  different  advance  ratios 
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Figure  4.14:  Velocity  profile  across  the  vortex  at  Azimuth  270 


Figure  4.15:  Flow  Behavior  Around  the  Blunt  Edge.  Qualitative  pictures  from  initial 
SPIV  data 
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Figure  4.16:  Streaklines  and  radial  velocity  contours  at  the  Blunt  Edge,  for  r/R=0.4, 
at  n  =  0.7  at  240  and  270  deg. 


Chapter  5 
Conclusions 


Below  we  integrate  the  findings  on  the  aerodynamics  of  rotor  blades  in  reversed  flow: 

1.  At  high  yaw,  a  fixed  wing  at  angle  of  attack  in  reverse  flow  develops  a  vortex 
on  the  sharp  edge. 

2.  The  variation  of  lift  and  pitching  moment  with  angle  of  attack  shows  substantial 
vortex  lift  effect. 

3.  The  lift  and  pitching  moment  variation  are  better-predicted  using  the  Polhamus 
leading  edge  suction  analogy  which  allows  accounting  of  the  vortex  flow  effect. 

4.  Surface  tuft  behavior  shows  vortex  bursting  near  the  blade  root  for  the  forward- 
yaw  reverse  flow  cases. 

5.  In  the  case  of  a  rotating  blade,  since  the  blade  is  usually  pitched  up  on  the 
retreating  side,  the  leeside  in  reverse  flow  is  the  bottom  side  of  the  sharp  edge. 
This  is  where  the  sharp  edge  vortex  forms. 

6.  At  240  degrees  azimuth,  a  strong  SEV  is  seen  to  grow  in  the  inboard  direction. 

7.  The  inboard  rise  in  vortex  strength  is  due  both  to  SEV  growth,  and  linear 
increase  in  onset  velocity. 

8.  Between  240  and  270  degrees,  the  SEV  bursts  and  detaches  from  the  surface  as 
the  yaw  decreases. 

9.  The  vortex  seen  at  270  degrees  azimuth  is  generated  earlier,  detached  and 
convects  with  the  blade. 
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10.  The  vortex  structure  at  300  deg  is  similar  to  that  at  270  deg. 

11.  Static  pressure  contours  in  the  flow  field  show  a  low-pressure  region  close  to 
the  surface  by  vortex  interaction,  resulting  in  vortex  lift. 

12.  Vortex  velocity  profiles  at  270  and  300  degrees  show  an  expanded,  solid-body- 
rotation  core  occupying  most  of  the  vortex  region  under  the  blade. 

13.  Surface  pressure  contours  at  270  and  beyond,  show  the  minimum  pressure 
occurring  0.15c  from  the  sharp  edge. 

14.  The  pressure  field  computed  from  the  velocity  field  explains  inward  radial  flow 
at  270  and  300  deg. 

15.  Radial  flow  is  outward  except  in  the  SEV  zone. 

16.  The  flow  around  the  blunt  edge  is  highly  3-D  and  mostly  attached.  At  240  deg. 
streamline  curvature  shows  negative  lift  generation  at r/R  =  0.5,  whereas  at 
270  deg.  and  r/R  =  0.5,  there  is  upward  mean  lift. 

17.  There  is  cycle-to-cycle  intermittency  in  the  blunt  edge  flow,  but  no  evidence  of 
periodic  separation. 
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